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ABSTRACT 



We have completed a 1.1 mm continuum survey of 7.5 deg 2 of the Perseus 
Molecular Cloud using Bolocam at the Caltech Submillimeter Observatory. This 
represents the largest millimeter or submillimeter continuum map of Perseus to 
date. Our map covers more than 30,000 31" (FWHM) resolution elements to a 
lcr RMS of 15 mJy/beam. We detect a total of 122 cores above a 5cr point source 
mass detection limit of O.18M , assuming a dust temperature of Tp = 10 K, 60 of 
which are new millimeter or submillimeter detections. The 1.1 mm mass function 
is consistent with a broken power law of slope ct\ = 1.3 (O.5M < M < 2.5M ) 
and a?2 = 2.6 (M > 2.5M ), similar to the local initial mass function slope 
Oi = 1.6 M < 1M , a 2 = 2.7 M > 1M ). No more than 5% of the total cloud 
mass is contained in discrete 1.1 mm cores, which account for a total mass of 
285M . We suggest an extinction threshold for millimeter cores of Ay ~ 5 mag, 
based on our calculation of the probability of finding a 1.1 mm core as a function 
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of Ay. Much of the cloud is devoid of compact millimeter emission; despite the 
significantly greater area covered compared to previous surveys, only 5 — 10 of 
the newly identified sources lie outside previously observed areas. The two-point 
correlation function confirms that dense cores in the cloud are highly structured, 
with significant clustering on scales as large as 2 x 10 5 AU. Our 1.1 mm emission 
survey reveals considerably denser, more compact material than maps in other 
column density tracers such as 13 CO and A v , although the general morphologies 
are roughly consistent. These 1.1 mm results, especially when combined with 
recently acquired c2d Spitzer Legacy data, will provide a census of dense cores 
and protostars in Perseus and improve our understanding of the earliest stages 
of star formation in molecular clouds. 

Subject headings: stars: formation — ISM: clouds — ISM: individual (Perseus) 



1. Introduction 

Observations of the earliest stages of molecular core collapse and protostellar formation 
are extremely important pieces in the puzzle of low mass star formation, not only illumi- 
nating the transition from dense cloud cores to later protostellar phases, but also providing 
clues about initial conditions and a crucial link between theoretical and empirical scenarios. 
Prestellar cores have no internal luminosity source and are therefore very cold (5 — 20 K), 
with spectral energy distributions (SEDs) that peak at submillimeter or millimeter wave- 
lengths. At very early times in an embryonic protostar's life, generally referred to as the 
Class phase (Andre et al. 1993), it is deeply embedded in an obscuring envelope of gas and 
dust. Most of the continuum emission from the hot young star is absorbed and re-radiated 
by the cool (10 — 30 K) dust envelope at far- infrared (far-IR) to millimeter wavelengths. 

We follow Gregersen & Evans (2000) and distinguish prestellar cores, which are expected 
to eventually form a star or stellar system, from starless cores, which are dense cores without 
IRAS sources that may or may not collapse in the future, by the presence of submillimeter 
or millimeter emission. Gregersen & Evans (2000) found that strong submillimeter emission 
correlates well with collapse signatures and the presence of NH3, making these sources likely 
to be prestellar. In this paper, we refer to a 1.1 mm core as any discrete source of 1.1 mm 
emission, a prestellar core as any 1.1 mm core lacking far-IR emission, and a protostar as 
any protostellar object with a substantial envelope (Class 0/1; see Lada (1987); Andre et al. 
(1993) for a description of classifications). 

The millimeter wavelength range is an ideal region in which to study such objects for 
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several reasons. At these long wavelengths thermal emission from dust becomes optically 
thin, and the observed flux density traces the total mass of dust in cores and envelopes. 
Because they lack an internal luminosity source, prestellar cores are generally invisible at 
shorter wavelengths. Additionally, flux density measurements in the millimeter contribute 
important information to the SED, helping to clarify the evolutionary state of a given object. 
Although deeply embedded objects have remained relatively elusive due to the difficulty of 
observing at submillimeter and millimeter wavelengths, continually improving technology 
and instrumentation are now making sensitive observations in this regime possible. 

Large format bolometer arrays like SCUBA (Holland et al. 1999), SHARC II (Dowell 
et al. 2003) and MAMBO (Kreysa et al. 1998) have made it feasible to scan relatively large 
fields at continuum submillimeter and millimeter wavelengths, while interferometers such as 
OVRO (Woody et al. 1998), BIMA (Welch et al. 1996), PdBI (Guilloteau et al. 1992) and 
the SMA (Ho et al. 2004) begin to resolve the details of individual sources. Nevertheless, 
time constraints have, for the most part, prohibited coverage of entire molecular clouds 
with current millimeter continuum instruments. Large, complete, high resolution surveys of 
entire star forming regions are necessary to minimize bias and systematic effects, improve 
statistics, and develop a clear framework within which to interpret the many observations of 
individual objects that are now available. Bolocam, a 144-element bolometer array designed 
for mapping large fields at millimeter wavelengths (Glenn et al. 2003), is well suited to the 
kind of large scale surveys now needed. 

Nearby molecular clouds such as Perseus, Ophiuchus and Serpens, where there is consid- 
erable evidence of ongoing star formation (e.g. Evans 1999), provide the best opportunity to 
observe stars in the earliest stages of their formation. Perseus is located in the larger Taurus- 
Auriga-Perseus dark cloud complex at about 3 A 30 m , +31°, and extending approximately 7° 
in Right Ascension (RA) and 3° in Declination (Dec). The Perseus cloud is often cited as an 
intermediate case between the low-mass, quiescent Taurus and turbulent, high-mass Orion 
star formation regions (e.g. Ladd et al. 1993, 1994), making it perhaps an ideal environment 
for studying "typical" low mass star formation. Two young clusters lie in Perseus: IC 348 
is a young infrared cluster of age about 2Myr containing several hundred members of total 
mass about 160 M (Luhman et al. 2003). NGC 1333 is a very young (< lMyr) highly 
obscured cluster with about 150 stars totaling ~ 79 M (Lada et al. 1996; Wilking et al. 
2004), and evidence of ongoing star formation. Perseus contains several Lynds and Barnard 
dark clouds, including B5 at the eastern end, Bl, and L1455 and L1448 at the western end. 
A number of dense ammonia (NH 3 ) cores have also been identified by Ladd et al. (1994). 

Recent estimates of the distance to Perseus range from 220 pc to 350 pc (e.g Cernis 
1990; Herbig & Jones 1983). Larger values (300 — 350 pc) are often adopted based on the 
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Perseus 0B2 association, which has a fairly well established distance of ~ 330 pc from 
Hipparcos parallax and reddening measurements of its members (e.g. Borgman & Blaauw 
1964; deZeeuw et al. 1999). There is evidence, however, that Per OB2 may lie well behind 
the complex of dark clouds in which we are interested, and are probably at a distance closer 
to 250 pc (Cernis 1993; Belikov et al. 2002). On the other hand, a single distance for the 
whole cloud might not be appropriate. Early CO observations (Sargent 1979) indicated a 
smoothly varying LSR velocity gradient across the cloud from v = +3km s _1 at the western 
end to v = +10km s~ x at the eastern end. Given this gradient in velocity, there may also be 
a distance gradient across the cloud. Extinction studies of several different regions point to 
an increase in distance from 220 pc to 260 pc moving from west to east (Cernis 1990, 1993; 
Cernis & Straizys 2003). 

Another possibility is the superposition of two or more clouds. Based on star counts, 
Cernicharo et al. (1985) suggest two dust layers at d ~ 200 pc and d ~ 300 pc. New data 
compiled by the COMPLETE team also indicate that Perseus may actually be composed 
of several separate clouds at different distances, projected together on the sky (Ridge et al. 
2005). We acknowledge that the structure and dynamics of the Perseus cloud are compli- 
cated; it is possible, perhaps likely, that it will ultimately be shown that Perseus is in fact 
a superposition of a number of smaller clouds. For the purpose of this work, however, we 
adopt a distance of 250 pc for the entire area surveyed, based on the most recent extinction 
studies and parallax measurements of IC 348 members (Cernis 1993; Cernis & Straizys 2003; 
Belikov et al. 2002). 

Perseus has been fully mapped in CO isotopes tracing densities up to a few thousand 
particles per cubic centimeter (Padoan et al. 1999; Ridge et al. 2005). Most previous submil- 
limeter and millimeter continuum and molecular line mapping of higher-density gas tracers, 
however, has been confined to the dense cluster region of NGC 1333 (e.g. Sandell & Knee 
2001), or focused on bright IRAS sources (e.g. Ladd et al. 1994) or energetic outflow sources 
(e.g. Froebrich et al. 2003; Mathews & Wilson 2002; Barsony et al. 1998), with the notable 
exception of the recent SCUBA survey by Hatchell et al. (2005). Conclusions about the 
cloud based on such existing observations may be problematic. The gas traced by CO obser- 
vations has relatively low density and therefore is not necessarily representative of current 
star formation activity, while isolated small-scale studies may be providing a biased picture 
of the cloud. 

The observations presented here have the distinct advantage that they overlap entirely 
with the Perseus fields observed with IRAC and MIPS as part of the Spitzer Legacy Project, 
"From Molecular Cloud Cores to Planet Forming Disks" ("Cores to Disks" or c2d; Evans 
et al. 2003). The c2d project has mapped five large star forming clouds with the IRAC 
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(A = 3.6 — 8.0 /tm) and MIPS (A = 24 — 160 fim) instruments on Spitzer. While millimeter 
observations are essential to trace core and envelope mass and directly observe prestellar 
cores, IR observations are necessary to characterize embedded protostars. Additionally, 
IR measurements are more sensitive to temperature than to total mass, making them an 
excellent complement to millimeter observations. 

We have completed millimeter continuum observations of the entire Perseus cloud with 
Bolocam at the Caltech Submillimeter Observatory (CSO). This map, observed at A = 
1.1mm during 2003, January- February, is the first unbiased, flux-limited survey of cores 
and protostars in Perseus at millimeter wavelengths. As part of the same project, the 
remaining two large c2d clouds accessible from Mauna Kea, Ophiuchus and Serpens, also 
have coordinated large scale Bolocam 1.1 mm and Spitzer observations (Young et al. 2005, 
in preparation; Glenn et al. 2005, in preparation). When joined with the Perseus data, this 
combined sample will provide a unique basis for comparing star forming properties in varying 
environments, without the systematic difficulties introduced by observations from different 
instruments or variable coverage. Detailed SEDs for this complete sample of objects will be 
made possible through the combination of Spitzer IR and Bolocam 1.1 mm fluxes, as well as 
850 /im fluxes when available, allowing the construction of a more quantitative evolutionary 
sequence than currently available, and calculation of statistical lifetimes. 

Here we present the results of our Bolocam survey of Perseus. In Section 2 we describe 
the observations, and in Section 3 we describe in detail the reduction techniques used for the 
Perseus Bolocam data, including an iterative mapping scheme we have developed to restore 
source brightness lost during sky subtraction. Results are presented in Section 4, including 
source flux, size, and mass statistics and a comparison of the cloud appearance at 1.1 mm to 
other column density tracers. In Section 5 we discuss the completeness limits of the survey 
and the mass versus size distribution, the 1.1 mm Perseus mass function, and the clustering 
characteristics of the cloud. We end with a summary and plans for future work in Section 6. 

2. Observations 

Continuum observations at 1.1 mm were made with Bolocam 1 at the Caltech Submil- 
limeter Observatory (CSO) on Mauna Kea, Hawaii during 2003, January 28-February 15. 
Bolocam is a 144-element bolometer array that operates at A = 1.1, 1.4, and 2.1 mm. The 
instrument consists of a monolithic wafer of silicon nitride micromesh AC-biased bolometers, 
cooled to 260 mK. During the Perseus observations 81 of 144 bolometers were operational. 



http: / / www.cso.caltech.edu/bolocam 
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The field of view (FOV) is T.h and the beams are well approximated by a Gaussian with 
a FWHM of 31" at A = 1.1mm. The focal plane is flat over the FOV due to a cold field 
lens coupling the telescope to the array (Glenn et al. 1998, 2003; Haig et al. 2004). All 
observations of Perseus were completed in the 1.1 mm mode, which has a bandwidth of 45 
GHz and band center at 268 GHz. 

The integrated intensity 13 CO map of Padoan et al. (1999) was used to define the area 
of the Bolocam observations, as shown in Figure 1. The chosen boundaries correspond to 
an approximate extinction limit of Ay ~ 2 mag (Evans et al. 2003), and were designed 
to overlap entirely with the area to be observed with Spitzer as part of the c2d Legacy 
project. Maps were made in raster scan mode in sets of three scans offset by —11", 0", 
+ 11", with a scan speed of 60 " s~ x and a subscan spacing of 162". A subscan is defined 
as one pass of the array across the field, and a scan as a set of consecutive subscans that 
cover the field entirely. Offset scans are necessary to obtain a fully sampled map because 
the array elements are separated on the sky by 1.5/A, whereas 0.5/A corresponds to Nyquist 
sampling. Simulations indicated that 11" offsets were optimal for obtaining the best coverage 
perpendicular to the scan direction. 

The total area was divided into three large rectangles for the most efficient scanning. On 
each of 19 nights Perseus was observed for approximately two hours before and after transit, 
with a total of 12 scans of each section completed. A number of subscans have been omitted 
from the final map due to bad weather, manifested as very large sky noise, or temporarily 
poor bolometer performance. The total observing time was about 40 hours for the 7.5 square 
degree region, with 31 hours actually spent integrating on-source, or an observing efficiency 
of 75%. This corresponds to a mapping speed for Perseus of 3 arcmin 2 mJy _2 hr _1 . No 
chopping was done for any of the observations, thereby retaining, in principle, sensitivity 
to large-scale structure up to the angular size of the array (7'5) . Chopping is not required 
because the bolometers are AC-biased, which elevates the signal band above the atmospheric 
and instrumental 1/f noise. Demodulation brings the signal down to near-DC, with the 
signal band determined by the scan speed and beam size. 

One night at the beginning of the run was devoted to pointing and calibration observa- 
tions only, in order to set the focus and pointing constants for the run. Each night secondary 
calibrators, including the bright Class object NGC 1333-IRAS4A, were observed approxi- 
mately every two hours, with calibrator sources from many areas of the sky used to derive a 
calibration curve for the entire run (see Section 3.3). All calibrator observations were taken 
at the same scanning speed as science fields (60" s _1 ), and most were small 4' x 4' scan maps. 
At least one primary calibrator, usually Mars, was observed nightly, and six large 10' x 10' 
beam maps made over the 3 week run. These beam maps are used to define the distortion 
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corrections and the beam shape, which is found to be very gaussian. Nightly sky dips were 
used to measure the sky and telescope optical loading (the quiescent optical power received 
by the bolometers from the sky, telescope, and dewar). Weather was mostly clear for the 
run, with an average zenith 225 GHz tau of r ~ 0.07, ranging between 0.05 and 0.09. 

3. Data Reduction 

Given that these observations utilize a new instrument, we describe the data reduction 
process for Bolocam data in general, as well as reduction techniques specific to the Perseus 
data, in some detail. Preliminary data reduction was accomplished using a reduction pipeline 
written by the Bolocam instrument team, as described below and in Laurent et al. (2005). 
A number of problems specific to the bright sources and very long subscans of more than 
a few degrees observed for this project required the development of additional reduction 
techniques, in particular the iterative mapping routine described in Section 3.4. 

Initial steps in the pipeline include calculating the pointing model as a function of az- 
imuth and elevation, and calculating the RA/Dec of each bolometer for every time sample. 
The bulk of the data reduction effort goes into removing the (considerable) sky noise. Sub- 
sequently, bright pixel spikes from cosmic rays are flagged, and the power spectral density 
(PSD) of each subscan is calculated. The N bo i Q bolometer timestreams are converted into 
a 2D pixel map using pixel offsets, subscan PSDs, and calibration constants. Finally, an 
iterative mapping scheme developed for this purpose is used to recover flux density lost in 
the sky subtraction process. 

3.1. Pointing 

All pointing observations made during the 2003, January observing run, including sec- 
ondary calibrators as well as additional variable sources such as quasars, are used to calcu- 
late precise pointing corrections. The resulting pointing model serves to refine the recorded 
telescope pointing position. The RMS of the global pointing model for 2003, January as 
determined from the positions of known galaxies in the Lockman Hole is 9.1" (Laurent et al. 
2005). The pointing model is somewhat better in Perseus, likely due at least in part to the 
proximity of many pointing sources to the Galactic plane. 

We compared our positions for known sources to the literature, using submillimeter, 
millimeter, or radio positions when possible, and otherwise IRAS positions. There was a 
small systematic offset in RA (<5RA = — 3'.'2 ± 1'.'3), but not large enough to warrant a 
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correction to positions, and no systematic offset in Dec (5Dec = — 0'/5 ± l'/5). We find a 
la RMS compared to previous positions of 7", independent of azimuth or elevation. Given 
that this dispersion includes potentially large uncertainties in the literature positions, as well 
as possible physical offsets between IR and 1.1 mm sources, we conservatively estimate the 
overall accuracy for point source positions in Perseus to be 7". In addition to uncertainties 
in source positions, pointing errors increase the effective beam size, causing sources to be 
blurred in the coadded map and affecting the measured size and peak flux density, but not 
the total flux density. 

Position offsets of the various bolometers in the array from the telescope pointing cen- 
ter, as well as optical distortions - collectively termed pixel offsets - are measured using 
fully sampled observations of planets. Distortion corrections account for distortion in the 
optics, due primarily to an off-axis ellipsoidal mirror, and secondarily to imperfect optical 
alignments. Corrections are typically of order 2" — 3". Pixel offsets of each bolometer are 
used to compute the RA and Dec value of every time sample in a subscan, and later used to 
convert bolometer timestreams into a 2D pixel map. 

3.2. Removal of Sky Noise 

The sensitivity of a given Bolocam observation is determined by the intrinsic sensitivity 
of and optical loading on the bolometers, the integration time, and the success of sky noise 
subtraction. The most important reduction step is the removal of sky noise, or cleaning. On 
scales comparable to or larger than the beam, the instrument is limited by 1/f noise, which 
is primarily atmospheric but also instrumental in nature. Sky noise originates as fluctuations 
in the brightness temperature, or column density, of the atmosphere and dominates most 
astronomical sources at A = 1.1mm. Because chopping is not done, this noise is present in 
the bolometer timestreams before cleaning. 

The bolometer beams overlap almost completely in the near-field where the sky noise 
originates, but do not overlap in the far-field where astronomical signals originate. There- 
fore, to first order the sky noise is identical for each bolometer and a sky template can be 
constructed quite simply by taking an average or median across the bolometer array at each 
point in time. This average cleaning method is appealingly simple, but does not deal well 
with multiple correlated 1/f noises with different correlation coefficients, or with 1/f noises 
that are correlated on spatial scales smaller than the array. The latter might arise if the 
beams do not overlap completely at the height of the sky noise. 

In principle, the correlated noise can be removed because it is correlated in time, whereas 
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the astronomical signal is correlated in space across the array. A more sophisticated approach 
that addresses this issue is Principal Component Analysis (PCA) cleaning (see Laurent et 
al. (2005) and references therein for a discussion of PCA cleaning). In a PCA analysis the 
raw timestreams are projected along eigenvectors, bringing out common modes, or principal 
components, in the data. Patterns common to all bolometers correspond to sky noise, so 
subtracting such common modes from the data is an efficient sky subtraction technique. Re- 
moving the first principal component is nearly equivalent to performing average subtraction. 

Any number of components can be subtracted from the data, each removing progres- 
sively less correlated 1/ / noise. The actual reduction in the RMS noise depends on the initial 
sky noise present, but typically removing 3 PCA components reduces the overall noise by 
10 — 30% compared to average cleaning. Although removing more components will reduce 
the noise further, the disadvantage of PCA cleaning is that higher components tend to re- 
move source flux density (most of which can be recovered, see Section 3.4). Tests performed 
on observations of Serpens with both compact and extended sources indicate that removing 
between 1 and 5 components is most effective at eliminating stripes from 1/f noise while 
retaining source flux density and structure (Glenn et al. 2005, in preparation). 



3.3. Mapping and Calibration 

To make a 2D pixel map from the bolometer timestreams, the pointing model and 
empirically derived pixel offsets are used to project each bolometer time sample onto an 
RA/Dec grid of pixel values. Timestreams are coadded using a weighted (by the inverse of 
the PSD) average. We bin the map at a resolution equal to 1/3 of the true instrumental 
resolution, or 10"/pixel, which gives sufficient hits per pixel for significant statistics without 
degrading the resolution. Given the nature of the instrument and observations, a single 
pixel in the map will contain data from many bolometers and many scans. We refer to the 
coverage map as an image of the number of hits per pixel, or seconds per pixel, in the map. 
The coverage is dependent on pixel size, scan strategy, the number of bolometers, and the 
number of scans in the map. The average coverage for Perseus is about 500 hits pix -1 (or 10 
s pix -1 ), varying by 30% across the map and leading to ~ 15% variation in the RMS noise. 

To maximize the signal to noise of, and thus our chance of detecting, point sources, we 
optimally filter the map. Because the signal from a point source lies in a limited frequency 
band, we can use an optimal (Wiener) filter to attenuate 1/f noise at low frequencies, as 
well as high frequency noise above the signal band. Attenuating the 1/f noise reduces the 
overall RMS per pixel by ~ y/3 (making the RMS/pixel ~ RMS/beam), maximizing the 
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probability of detection for faint point-like sources. The optimal filter g(q) is given by: 

SW)\ 2 IJ(q)<Pq 

where q is the spatial wavenumber, J(q) is the azimuthally averaged PSD, s(q) is the Fourier 
transform of the beam, and g(q) is normalized so that the peak brightness of point sources 
is preserved. Note that the optimally filtered map is used for source detection only; all maps 
shown are unfiltered. 

The flux calibration factor at any given time depends on the bolometer DC resistance, 
which changes with atmospheric attenuation and bolometer optical loading. Therefore, the 
calibration cannot be applied as a single factor to the final map, but must be calculated as a 
function of time for each subscan based on the average DC resistance of that subscan. This 
is a powerful calibration method because it makes real-time corrections for the atmospheric 
attenuation and bolometer operating point using the bolometer optical loading, effectively 
providing flux calibration on the timescale of minutes. 

In practice, the calibration (in mV/Jy) as a function of the DC voltage is determined us- 
ing a second-order polynomial fit to observations of Mars, Saturn, and a number of secondary 
calibrator sources (see Laurent et al. (2005) for an expanded discussion). All calibrators with 
reliable known flux densities (one observed every 1-2 hours per night throughout the run) are 
used to define the calibration curve, not only those observed during a given science obser- 
vation. This enables a comprehensive calibration curve over a range of atmospheric optical 
depths. Reference flux densities are from the JCMT planetary flux website 2 and Sandell 
(1994), corrected for the Bolocam bandpass. The effects of non-linearity due to optical load- 
ing from the sky and the finite beam size are also accounted for, using large beam maps of 
Mars and assuming a uniform disk model. Measured flux densities are from a gaussian fit 
after application of a point source filter. 

The resulting peak-normalized map can be divided by the beam area to obtain a surface 
brightness-normalized map. The integral of a point source in the surface brightness map 
returns the total flux of the point source, thus this map is suitable for photometry. The 
absolute calibration uncertainty, derived from the deviation of measured calibrator values 
from the quadratic fit, is 9.7% (la; Laurent et al. 2005). 



2 http://jach/hawaii.edu/JACpublic/JCMT/software/bin/fluxes.pl 
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3.4. Iterative Mapping 

3.4.I. Method 

Although they utilize different methods, both average and PCA cleaning contain a step 
that essentially removes the mean of each bolometer subscan (see Section 3.2). This step 
is necessary to eliminate sky noise, but when there is a bright source in a subscan it biases 
the mean. Consequently, sky subtraction introduces negative lobes around bright sources, 
which are asymmetric in the scan direction. Furthermore, subscans containing sources tend 
to be under-weighted in the coadded map because the source brightness contaminates the 
integrated PSD, causing a decrease in the weighting factor. Both effects tend to suppress 
source flux density and are mildly dependent on the brightness and structure of the source. 
There is an additional effect due to PCA cleaning that non-linearly removes the flux density 
of bright sources as more PCA components are removed. To correct for diminished source 
brightness and negative artifacts introduced by the above effects, we have developed an 
iterative mapping code that robustly restores lost flux density and structure to the map. 

The iterative mapping algorithm we have implemented iteratively subtracts a source 
model from the real data (somewhat similar to CLEAN (Hogbom 1974; Schwarz 1978), but 
working in the image plane). The following is a more robust method than using, for example, 
a source model comprised of the sum of many gaussians, because many of our sources are 
likely to be extended and non-gaussian. For the following, j refers to bolometer number 
(j = 1 — 81), and i refers to the iteration number (i = — N, where the zeroth iteration 
% = indicates raw or cleaned data before any source model has been subtracted). We begin 
with the raw timestream data, t/,i=o- These data are sent through the cleaning and mapping 
process to produce the zeroth iteration map, M i=Q , which contains negative artifacts and is 
missing some fraction of the flux density of each source. Next, a cut is applied at +N a a to 
M i=0 , removing any negative pixels as well as most of the noise (pixels with values < +N (T 
times the RMS noise are set to zero). We now have a source model map, M[, our current 
best guess of the true source flux density. 

From the source model M[ a model timestream is generated for each bolometer 
and subtracted from the raw timestreams, dtj ti = tj >0 — t'-y The difference timestreams 
dtj t i contain residual source flux density that was missing from the original map. When the 
difference timestreams are subsequently cleaned and mapped to produce a residual map dM iy 
there is much less contamination of the sky template by source brightness, so the negative 
artifacts are greatly reduced. A threshold cut at +N a a is again applied to dMi producing 
a residual source model dM^, which is added to the original source model to create a new 
source model for the next iteration {M' i+1 — M[ + dM[). This process is iterated until there 
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is no remaining residual source flux density in the difference timestream. After N iterations, 
the final residual map dM^, containing only noise and any source flux below the threshold 
cutoff, is added to the last source map to create the final map (M N = + dM N ). 

3.4-2. Performance 

The performance of iterative mapping depends on two parameters: the number of PCA 
components removed during cleaning (Npca), and the RMS threshold cut (N a ) used to make 
the source model. While iterative mapping greatly improves the appearance and photometric 
reliability of the map, our chosen values of Npca and N c will introduce some systematic 
uncertainties into the final photometry. Using Monte Carlo simulations, which consisted of 
inserting simulated sources into the real map before the iterative mapping process, we found 
that Npca = 3 and N c = 2 produce the best combination of accurate final photometry 
(recovery of lost source flux density) and reduced noise compared to average cleaning (by 
~ 20 — 30%). Using a value oi N a < 2 tends to introduce noise into the source model, 
increasing the RMS noise and raising the background level of the map; increasing N a leaves 
source flux density out of the source model, leaving some negative artifacts and making it 
difficult to recover extended source features. Decreasing N P qa leaves significantly more 1/ / 
noise in the map and increases the overall RMS noise; if we increase Npca much above 3 it 
again becomes difficult to fully correct for the loss of source brightness, which is severe and 
non-linear for Npca > 5. 

Figure 2 demonstrates the results of iterative mapping on a portion of the Bolocam 
map containing the crowded region NGC 1333. The left panel shows NGC 1333 after PCA 
cleaning only with Npca — 3, where the dark blue areas are negative bowls caused by 
cleaning. The same region after 10 iterations is shown on the same intensity scale in the 
center panel. All sources have increased considerably in brightness (the brightest peak by 
14%), and most of the negative features have been corrected. Although some negative 
artifacts remain, their extent and intensity is greatly reduced; the worst remaining negative 
pixel in the iterated map is —88 mJy/beam, compared to —238 mJy/beam in the cleaned 
map. For comparison, the 850 /xm map of the central portion of NGC 1333 from Sandell & 
Knee (2001) is also shown on approximately the same intensity scale. Clearly, the extended 
structure in the higher resolution (14") 850 /zm map is recovered in the 10 iteration Bolocam 
map. The recovered structure is almost certainly real, therefore, and we are not missing 
any structure at 1.1 mm that is present at 850 /iin (although Sandell & Knee (2001) note 
that flux densities for extended emission in their map are unreliable due to similar residual 
artifacts). 
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Monte Carlo tests were done to quantify the effectiveness of the iterative mapping and 
characterize any systematic errors introduced throughout the cleaning and iterative mapping 
process, using N PC a = 3 and N a = 2. Simulated gaussian sources of varying amplitudes 
(3 — 300(x) and sizes (30" — 250" FWHM) were inserted into an empty region of the map, 
then cleaned and iterated. Finally, the resulting peak amplitude, flux density and FWHM 
were measured. The fractional missing peak flux density after 10 iterations for a range of 
test sources is shown as a function of the input peak amplitude in Figure 3. Sources of 
FWHM 30", 50", 80" and 120" are indicated. The la RMS noise per beam is also plotted as 
a percentage of the input amplitude. The fact that most points lie beneath this line indicates 
that the residual bias in measured flux densities after iterative mapping are consistent with 
the RMS noise. Only for the largest sources (120" FWHM, few of which are found in the 
real map), is there a residual systematic reduction in flux density that is larger than the 
RMS. The measured FWHM, axis ratio, and position angle are not significantly affected by 
either cleaning or iterative mapping. 

Figure 4 shows the fractional lost flux density in a 40" aperture as a function of iteration 
number for four representative FWHM sizes. A range of input amplitudes from S/N = 5 
(detection limit) to S/N = 175 are plotted as different line styles. Large, bright sources 
(bottom right panel) are the most affected by cleaning - more than 50% of the flux density 
is removed during PCA cleaning - but are almost entirely corrected by iterative mapping, 
with only ~ 2% residual missing flux density after 5 iterations. For large faint sources < 20<7 
(e.g. the dotted line, bottom right), iterative mapping is unable to recover all the flux density 
to better than 5 — 10%, but note that the la RMS in a 40" aperture is > 10% for sources 
< 18a. Sources larger than ~ 200" FWHM (not shown) are not fully recovered even after 
20 iterations, making this the effective maximum source size detectable. No sources larger 
than 150" FWHM are measured in the real map. 

We conservatively assign a systematic uncertainty of 5% to all integrated flux densities. 
In almost all cases the measured flux density is lower than the true value, making the 
uncertainty a bias in that measured values underestimate the true flux densities. We do 
not attempt to correct for this residual bias, instead including the 5% iterative mapping 
uncertainty in the overall systematic uncertainty of 15%. 

4. Results 

The final coadded 10"/pixel Perseus map is shown in Figure 5. The map has a beam 
size of 31" and a total area of 7.5 square degrees, covering a total of 3.4 x 10 4 resolution 
elements. The mean RMS is 15 mJy/beam. Well known regions are marked, including the 
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conspicuous bright cluster NGC 1333. Note that the infrared cluster IC 348 actually lies 
slightly to the northeast of the group of 1.1 mm sources indicated as IC 348 here. 

While nearly all previously known sources are seen, with the exception of some crowded 
regions where blending occurs (e.g. in NGC 1333 and L1448), perhaps the most striking 
feature of the map is its relative emptiness. In fact, few new sources are detected far from 
known cluster regions, and those that are tend to be faint objects near the detection limit. 
There does not appear to be any symmetric extended emission > 3' in the map. Although 
it should in principle be possible to recover symmetric structures up to the array size (7f5), 
our simulations show that sources > 4' are severely affected by cleaning and difficult to fully 
recover with iterative mapping. 

4.1. Source Identification 

Discrete sources of 1.1 mm emission, or cores, are identified within the optimally filtered, 
iterated map using a simple peak finding algorithm. An optimal filter (as described in 
Section 3.3) is applied to the iterated map to facilitate the detection of point sources. This 
decreases the noise in the map, as well as increasing the peak brightness (and therefore the 
probability of detection) of sources larger than the beam. After optimal filtering the noisy 
edges of the map are removed using a cutoff based on the coverage map. A cutoff at 20% of 
the peak coverage is chosen empirically based on the number of false edge sources detected. 
All peaks 5a above the average RMS are flagged as possible sources. 

To make the final cut, candidate sources must also be 5a above the local RMS, and 
have a well defined centroid. The centroid is a weighted average position based on the 
surface brightness within a specified aperture. The local RMS noise per beam is calculated 
in small (~ 45 arcmin 2 ) boxes in a noise map, from which sources have been removed using 
the source model generated during iterative mapping (see Section 3.4). The average RMS 
is 15 mJy/beam, varying across the map by 15%. Most of the variation in RMS is due to 
30% variations in observational coverage rather than to a change in the calibration, which is 
very stable. Peaks separated by more than one beam size from a previously identified source 
centroid are considered to be distinct sources. A few false sources at the edge of the map 
were not automatically cut and had to be removed by visual inspection. 

A total of 122 confirmed sources are identified in the optimally filtered map, the locations 
of which are indicated by small circles on the unfiltered map in Figure 6. This figure also 
includes insets of the densest source regions. Because the RMS varies very little across the 
map, the lack of sources over large regions of the image is real. Many of the sources seen in 



- 15 - 



the 1.1 mm Bolocam map were very recently identified at 850 /im in the SCUBA survey by 
Hatchell et al. (2005). The number of new sources not previously identified at submillimeter 
or millimeter continuum wavelengths is ~ 60/122, most of which are within the region 
covered by Hatchell et al. (2005) but not detected in that survey due to their somewhat 
higher mass detection limit (0.4 M© (12K) vs. 0.2 M & (10 K) for the present survey). A 
number of these new sources lie within dense NH3 cores (Per3-Per9; Ladd et al. 1994), or 
coincide with IRAS sources. Most are faint, and in the vicinity of known groups, although 
there are a few more isolated sources. A selection of new sources chosen to demonstrate the 
wide range in source properties is shown in Figure 7. 

We additionally identify cores using the 2D Clumpfind routine of Williams et al. (1994). 
Clumpfind is useful for separating sources in crowded regions, where it may be more effective 
than aperture photometry in defining total flux densities. It assigns pixels to each source 
by first contouring the map with a small contour interval (2a). For each peak, contours are 
followed using a "friends-of-friends" algorithm down to one contour below the 5a detection 
limit, or until they run into another source at the same contour level. Clumpfind makes no 
assumptions about the clump profile, so the sources need not be gaussian. 

Clumpfind identifies 119 cores in the unfiltered map. Hereafter we refer to the sample 
of cores detected using peak finding in the optimally filtered map as "Peak-find sources", 
and the sample identified with Clumpfind as "Clumpfind sources". The total number of 
Clumpfind and Peak-find sources is very similar, but there is not a one-to-one correlation 
between the two sets. Fewer faint sources are found with Clumpfind because sources are 
detected above 5a using the RMS per pixel of the unfiltered map, which is ~ times the 
RMS per beam used to detect sources in the optimally filtered map. On the other hand, 
Clumpfind breaks up many of the brighter Peak-find sources into multiple sources. Faint 
extended sources that one would consider single if examining by eye are often partitioned 
into multiple sources by Clumpfind, and for this reason we favor the Peak-find source identi- 
fication and photometry. We keep Clumpfind calculations for comparison to previous work, 
as well as to illustrate differences between the two methods. 



4.2. Comparison to Molecular & A v Maps 

Large-scale CO and extinction maps of Perseus are available for public use as part 
of "The COMPLETE Survey of Star-Forming Regions" 3 (COMPLETE; Goodman et al. 
2004). These observations have been coordinated to cover the Spitzer c2d area of Perseus, 



3 http://cfa-www.harvard.edu/COMPLETE/ 
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therefore they also overlap with the Bolocam observations presented here. For the following 
we make use of the COMPLETE 13 CO map, and calculate our own NIR extinction (Ay) 
map. It is most instructive to compare the 1.1 mm, Ay, and 13 CO maps if they are all 
converted to the same resolution and column density scale. The three maps are shown as Ay 
contours overlaid on the grayscale 1.1 mm emission in Figures 8-10. Contours in all plots are 
Ay = 2, 3, 4, 6, 8 mag except for the 1.1 mm data (Figure 10), where Ay = 2, 5, 10, 15, 20 mag 
to avoid confusion, and smoothed to an effective resolution of 5'. 

The extinction contours in Figure 8 were calculated from the H-K s color excesses of 
2MASS sources using the NICE method (e.g. Lada et al. 1999; Huard et al. 2006) and 
convolving with a Gaussian beam with FWHM of 5'. This method depends on background 
stars to probe the column density through the cloud. In constructing the extinction map, 
we eliminate from the 2MASS catalog most foreground and embedded sources that would 
yield unreliable extinction estimates. Foreground sources were identified by color excesses 
representing a small line-of-sight visual extinction, Ay^os < 3 mag, within a 5' x 5' region 
exhibiting a large mean extinction of Ay > 8 mag. Embedded sources were identified as 
those sources positioned at least la redward of the reddened main-sequence region on a (J-H, 
H-K) color-color diagram. We use only those sources with reasonably good photometry at 
J, H, and K s (photometric errors less than 0.5, 0.2, and 0.2 mag, respectively) to construct 
the extinction map. The color-excesses were computed assuming that the intrinsic color of 
non-extincted stars in off-cloud fields are representative of the background star population. 
These two off-cloud regions, having a total of more than 5700 stars, were 1?5 x 0?4 fields 
centered on a = 03 /l 56 m 48.0 s , 5 = 35°06'00" and a = 03 fe 49 m 12.0 s , 5 = 27°54'00" (J2000.0). 

The NIR-derived extinction contours in Figure 8 verify that the Bolocam 1.1 mm map 
covers most of the Ay > 2 mag area. All bright 1.1 mm sources occur in regions of high 
extinction, but the reverse is not necessarily true; high extinction does not guarantee the 
presence either of young protostars or dense millimeter cores with strong 1.1 mm emission. 
Note that the extinction map saturates around Ay > 10 mag. 

Where 13 CO is optically thin it provides another measure of column density. This may 
only be true for a small range of densities, however, as 13 CO can be under-excited at small 
densities and becomes optically thick at Ay ~ 6 mag for a typical linewidth of Av ~ 1 km s -1 . 
To convert 13 CO integrated intensity to H2 column density, we assume LTE, r< 1, and an 
excitation temperature T x = 10 K, following Dickman (1978). Ay contours are calculated 
taking A" iTB ( 13 CO) = 2.5 x 10 15 (A v -0.8) cm~ 2 in Perseus for 1 < A v < 5 mag (Bachiller & 
Cernicharo 1986), and shown in Figure 9. The Ay( 13 CO) and NICE extinction maps display 
the same general morphology. Although the 13 CO contours do not reach peak densities as 
high as those in the NICE Ay map, they do become more compact compared to the NIR 
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extinction derived Ay, especially around 1.1 mm sources. 

H 2 column density is calculated from the 1.1 mm map using a conversion from 1.1 mm 
flux density per beam (S^™ m ) °f 

obeam 

# (H 2 ) = = v - (2) 

Here Vl beam is the beam solid angle, m H is the mass of hydrogen, K Llmm = 0.0114 cm 2 g _1 is 
the dust opacity per gram of gas, B u is the Planck function, To=10 K is the dust temperature, 
and hh 2 — 2.8 is the mean molecular weight per H2 molecule, which is the relevant quality for 
conversion to extinction. We extrapolate K\,\ mm from Ossenkopf & Henning (1994, Table 1 
column 5, hereafter OH5) for dust grains with thin ice mantles, coagulated for 10 5 years at 
a gas density of 10 6 cm" 3 . A gas to dust ratio of 100 is assumed. OH5 dust has found to be 
the best choice for star-forming cores by several authors using radiative transfer modeling 
(Evans et al. 2001; Shirley et al. 2002; Young et al. 2003). Column densities calculated 
from thermal emission depend on To, which varies with position. Thus there may be local 
discrepancies with other tracers when T D is not independently known, but agreement should 
be good overall. 

We assume a conversion from column density to Ay of 7V(H 2 )/Ay = 0.94 x 10 21 
cm 2 mag -1 (Frerking et al. 1982). The resulting A V {1.1 mm) contours in Figure 10 clearly 
demonstrate that the 1.1 mm emission is considerably more compact compared to the other 
column density tracers. It appears that only the densest material is traced by the 1.1 mm 
emission, as manifested in both the compact nature of the Ay (1.1 mm) contours and the 
high Ay reached in the bright cores {Ay^ max ~ 17 at 5' resolution). The 5<r detection limit 
per 31" beam corresponds to a minimum Ay ~ 5 mag, so it is not surprising that there 
is not 1.1 mm emission seen in the lower Ay regions of the Ay and 13 CO maps. We note, 
however, that our simulations indicate that the Bolocam map may not be sensitive to more 
extended emission (> 4') even if it is present in the cloud. 



4.3. Source statistics 

Source identifications, positions, peak flux densities (S pea k), and signal to noise (S/N) 
for the 122 1.1 mm sources identified in the Bolocam map are listed in Table 1. The S/N is 
measured in the optimally filtered map, because this is the S/N that determines detection. 
Photometry and all other source properties are measured in the unfiltered, surface brightness- 
normalized map. Sp ea k is the peak pixel value in mJy/beam as measured in the 10"/pixel 
map. The uncertainty in S pea k is the local RMS/beam, calculated as described in Section 4.1. 
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There is an additional systematic error of ~ 15% associated with all flux densities, from the 
absolute calibration uncertainty (10%) and the systematic bias remaining after iterative 
mapping (~ 5%, see Section 3.4). Table 1 also lists the most commonly used name from the 
literature for known sources, and indicates if the 1.1 mm source is coincident (within 40") 
with a MIPS 24 /im source from the c2d database (Rebull et al. 2005, in preparation). The 
c2d MIPS data for some regions of Perseus are not yet available, in which case "N/A" is 
listed in the table. For these sources, the presence of an IRAS source (within 1') is indicated 
if appropriate. 

Table 2 lists photometry in fixed apertures of diameter 40", 80", and 120", the total 
integrated flux density (S v ), total mass, peak Ay, full- width at half maximum (FWHM) 
sizes along the major and minor axes, position angle (PA, measured east of north), mean 
particle density (n), and morphology descriptions for each source. Integrated flux densities 
are measured assuming a sky value of zero, and corrected for the gaussian beam so that 
a point source has the same integrated flux density in all apertures with diameter greater 
than the beam. No integrated flux density is given if the distance to the nearest neighboring 
source is smaller than the aperture diameter. 

To calculate the total flux density, we compute integrated flux densities in apertures of 
30" — 160" in intervals of 10", and use the largest aperture diameter that is smaller than the 
distance to the nearest neighbor. An aperture of 120" is generally sufficient unless the source 
FWHM is > 100", in which case apertures up to 160" are used. If the flux density decreases in 
larger apertures (due to residual negative artifacts), then the aperture giving the largest flux 
density is used. Uncertainties in all integrated flux densities are a ap = (? beam ^/ A ap /A beam , 
where a beam is the local RMS/beam and (A ap , A beam ) are the aperture and beam areas 
respectively, and do not include an additional 15% systematic uncertainty. The FWHM and 
PA are from the best fit elliptical gaussian after masking out nearby sources using a mask 
radius equal to half the distance to the nearest neighbor. The errors given are the formal 
fitting errors and do not include uncertainties due to noise or residual cleaning effects, which 
are of order 10 - 15% for the FWHM and ~ 5° for the PA. 

The distribution of source flux densities is shown in Figure 11. In addition to peak 
flux densities (Jy/beam) and total flux densities (Jy) from aperture photometry of Peak- 
find sources, Clumpfind 3<r flux densities (Jy) are also shown. Clumpfind flux densities are 
calculated from the source map generated by Clumpfind as the sum over all > 3<x pixels 
assigned to a given source, divided by the beam area. We find empirically that Clumpfind 
is good at separating, and calculating the correct flux density for, bright, crowded sources. 
It is not, however, very effective at determining the total flux density of isolated faint or 
extended sources. For example, faint sources will have a smaller total flux density than for 
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aperture photometry because Clumpfind only includes pixels above 3a. 

The large average size of sources accounts for the fact that the mean total flux density 
(0.96 ± 1.21 Jy; hereafter, numbers quoted are mean ± dispersion in the sample, not mean 
± error in the mean) is significantly larger than the mean peak flux density (0.35 ± 0.56 
Jy/beam). Clumpfind sources tend to have slightly larger 3a flux densities (1.13 ± 1.80 Jy), 
in part because bright, crowded sources may be integrated over a larger area than with 
aperture photometry, and in part because one Clumpfind source sometimes encompasses 
more than one Peak-find source. 

Figure 12 gives the distribution of major and minor FWHM sizes, as well as the full- 
width at 3a (FW3<r) sizes of Clumpfind sources, defined as FW3<r = 2 a/ N P i X A pix /7r, where 
N pix is the total number of pixels assigned to a source and A pix is the pixel area. The average 
minor axis FWHM is 58" ± 17", the average major axis FWHM is 80" ±27", and the average 
mean FWHM is 68" ± 20". Clumpfind sources have a large average FW3a of 77" ± 30", as 
expected given that the size is measured at 3a rather than at half-max. There are very few 
sources near the resolution limit of 31". Most sources are somewhat elongated as well as 
extended; the average axis ratio is 1.4, and the distribution (Figure 13) extends to large axis 
ratios > 2. Note, however, that measured axis ratios < 1.2 are found to be unreliable based 
on simulations, and should be considered indistinguishable from an axis ratio of unity. 

The total mass M of gas and dust in a core is proportional to the total flux density S u , 
assuming the dust emission at 1.1 mm is optically thin and both the dust temperature and 
opacity are independent of position within a core: 

d 2c l 

M — \ , (3) 

B v {TdW 

where K\ \ mm = 0.0114 cm 2 g" 1 is the dust opacity, d = 250 pc is the distance, and To is the 
dust temperature. Although the millimeter emission arises only from the dust, we can infer 
the total mass of gas and dust by assuming a gas to dust mass ratio of 100, which is included 
in /ti.imm- For the masses in Table 2, we assume a single dust temperature To = 10 K for 
all sources. The uncertainties given are from the uncertainty in the total flux density only. 
Other sources of error from k, T d , and d (up to a factor of 4 or more) are discussed in 
Section 5.2. 

For dense regions without internal heating, the mean temperature is about 10 K, warmer 
on the outside and colder on the inside (Evans et al. 2001). Centrally heated cores will be 
warmer on the inside, but much of the mass is located at low temperatures. Shirley et al. 
(2002) and (Young et al. 2003) found good agreement with masses of Class and Class I 
sources determined from detailed models using To ~ 15 K. Taking 10 K is a reasonable 
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compromise to cover both prestellar and protostellar sources, but keep in mind that it will 
overestimate the masses of the latter by a factor of 2 — 3. 

The peak A v in Table 2 is calculated from the peak 1.1 mm flux density S peak as in 
Equation 2. The average peak A v of the sample is 24.7 mag. The mean particle density 
for each source is estimated as (n) = M/(4/3irR 3 fi), where M is the total mass, R is the 
mean deconvolved HWHM size, and n P = 2.33 is the mean molecular weight per particle. 
The average (n) of the sample is 4.3 x 10 5 cm -3 . The morphology keywords listed indicate 
if the source is multiple (within 3' of another source), extended (major axis FW at 2a > V), 
elongated (axis ratio at 4<r > 1.2), round (axis ratio at 4a < 1.2), or weak (peak flux 
densities less than 5 times the RMS per pixel in the unfiltered map). Monte Carlo simulations 
(Section 3.4) indicate that we cannot recover structures larger than ~ 200", and we do not 
find any sources larger than ~ 120" in the real map. We do not resolve any source pairs closer 
than 32", close to the minimum separation of 30" required by the peak finding algorithm. 



5. Discussion 

5.1. Completeness & the mass vs. size distribution 

The distribution of total mass vs. FWHM size is shown in Figure 14. The solid line 
shows the expected mass detection limit for gaussian sources assuming a simple scaling with 
source area. The gaussian mass limit varies with source size {Mu m oc size 2 ) because our 
ability to detect sources is based on the peak flux density above the noise (5a), whereas 
the mass is calculated from the total flux density (S u oc size 2 ). The real mass completeness 
limit is more complicated, even for gaussian sources, due to our reduction and detection 
techniques. We also show, therefore, empirical mass detection limits as a function of size for 
10%, 50% and 90% completeness. 

Completeness is determined by Monte-Carlo simulations, taking into account effects 
from cleaning, iterative mapping, and optimal filtering. Simulated sources of varying total 
mass and size are inserted into an empty region of the real map before cleaning and iterative 
mapping. The 10% completeness limit is the mass at which 10% of the simulated sources 
are detected in the optimally filtered map. Nearly all detected sources lie above the 10% 
completeness limit, as expected. Note that the completeness curve represent the true, not 
measured, mass and size of input simulated sources. 

Typical measurement error bars in M and FWHM are shown for 50" and 100" FWHM 
sources near the detection limit. For these, <tm are from the uncertainty in the integrated 
flux, and ctfwhm are estimated from simulations. The maximum size of the pointing-smeared 
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beam is indicated by a shaded band. Using the optimally filtered map to detect sources 
tends to lower the mass limit for sources with FWHM > 30" because the peak is enhanced 
by the filter. Conversely, very large sources (> 100") tend to have a higher mass limit 
because they are not fully recovered by iterative mapping. Both effects are illustrated by 
the empirical 50% completeness curve, which falls below the gaussian 50% completeness line 
for FWHM ~ 40" - 60", and above for FWHM > 100". 

There is an additional complication due to the fact that the sources may not be gaussian 
in shape. Without knowing the true source structure this effect is difficult to quantify. We 
explore the effects of different source structures by running completeness tests for sources 
with a Bonnor-Ebert (BE) sphere profile (Ebert 1955; Bonnor 1956), which has been found 
by several authors to be a good representation of the structure of some prestellar cores 
(Schnee & Goodman 2005; Kirk et al. 2005; Evans et al. 2001; Ward-Thompson et al. 1994). 
The 50% completeness limits for BE spheres are indicated in Figure 14 by a thick shaded 
line. FWHM sizes of BE spheres, measured with a gaussian fit for consistency with other 
limits, are set by the chosen outer radius (r Q ) in the BE model. The total mass is determined 
by a combination of r Q and the central density n c , and detection primarily by n c . The BE 
sphere models used to calculate completeness have n c = (8 x 10 4 ,8 x 10 4 ,9.5 x 10 4 ) cm~ 3 , 
r = (8 x 10 3 , 1.5 x 10 4 , 3 x 10 4 ) AU, and M tot = (0.45, 0.91, 1.89) M Q . 

Figure 15 is similar to Figure 14 except that the mass versus FW3a size is plotted 
for both Peak-find and Clumpfind sources. For Peak-find sources, FW3a is scaled from the 
FWHM assuming a gaussian profile (FW3a =FWHMx ^ln (S , peafc /3<T)/ln2). Masses for the 
two are calculated slightly differently: the Peak-find mass is the total mass from aperture 
photometry, whereas the Clumpfind 3<r mass is from the integrated flux density within the 
3o- contour. The total mass for Peak-find sources (circles) has not changed from Figure 14, 
but the distribution of mass vs. size looks very different because the FW3a size tends to be 
larger than the FWHM for bright sources (by about a factor of 2, because the 3a contour 
is well below the half-max), but smaller than the FWHM for faint sources. The apparent 
decrease in scatter is not real, but simply an effect of the size definition used. For faint 
sources the Clumpfind 3<r mass is smaller than the Peak-find aperture photometry mass, 
because Clumpfind only integrates down to 3<x, which corresponds to a small aperture for 
faint sources near the 5a detection limit. 

The appearance of the mass vs. size distribution depends very strongly on the definitions 
of both size and total mass, as illustrated by Figures 14 and 15, making comparisons between 
plots calculated in different ways quite deceptive. For example, the 3a mass vs. 3a size 
distribution for Clumpfind sources in Figure 15 (crosses) is very similar to that in Ophiuchus 
from Johnstone et al. (2000), but comparing the Ophiuchus plot to the Peak- find total mass 
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vs. FWHM size would make the two distributions seem very different. 

We find it significant that there seem to be very few point sources in Perseus at 1.1 mm. 
This is most clearly illustrated in Figure 14, which demonstrates an obvious paucity of both 
faint and bright sources between FWHM~ 30" — 50" (although there are low mass, compact 
sources in Figure 15, these are artificially created by the 3<r cutoff). Given that the pointing- 
smeared beam is no larger than 34" (for pointing errors < 7") and the average FWHM is more 
than twice the beam size (68" vs. 31"), the majority of sources are significantly extended. In 
fact, the average deconvolved size is 61" = 1.5 x 10 4 AU. This mean size is inconsistent with 
descriptions of cores as truncated spheres, which have been used to model very low mass 
cores. A truncated power law model for an 0.3 M core requires an outer radius of 2 x 10 3 AU 
(Young & Evans 2005), which would be a point source at the resolution of Bolocam. Our 
mass limit is 0.18 M for a point source at T D = 10 K, so we should be sensitive to such 
compact, low mass cores if they are present in Perseus. 



5.2. The 1.1 mm Mass Function 

The differential mass function dN/dM for all 122 1.1 mm sources is shown in Figure 16. 
The average mass of the sample is (M) = 2.3M , with individual masses ranging from 0.2 
to 26M Q . A dust opacity of K, 1Amm = 0.0114 cm 2 g" 1 (OH5), temperature of T D = 10 K, 
and distance of d = 250 pc are assumed for all masses. Error bars are \fN statistical 
errors, and are meant to demonstrate the typical uncertainties from photometry, but do not 
include errors due to uncertainties in the distance or dust properties. Completeness becomes 
important around 0.2 M & (50% completeness) for point sources, and around 0.8 M & for 
sources with a FWHM of ~ 70", the average FWHM of the sample. The mass distrubution 
is not corrected for incompleteness, therefore the turnover below ~ 0.8 M may not be real. 

The best fit to the mass function is a broken power law 

1 * oc M~ a (4) 



dM 



with cti = 1.3 ± 0.3 (O.5M < M < 2.5M ) and a 2 = 2.6 ± 0.3 (M > 2.5 M ). This result 
has a reduced chi-squared of x 2 = 0.4 and is true for any break mass between 2M and 
3M . The best fit broken power law is shown in Figure 16. The best fit single power law 
(a = 2.1 ± 0.1, M > 0.8 Mq) is also shown, although this is not as good a fit to the data 
(X 2 = 0.8) 

Johnstone et al. (2000) report a x = 1.5 (M < 1 M ) and a 2 = 2.5 (M > 1 M ) for 
850 /jm sources in the central 700 arcmin 2 of Ophiuchus assuming a single dust temperature of 
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20 K and using Clumpfind photometry. The Oph slopes are very close to those for the Perseus 
sample, and if we assume the same dust temperature as Johnstone et al. (2000) (T = 20 K), 
the break mass of ~ 2.5 M becomes ~ 1 M (the shape of the mass function does not change 
with Td as long as a single temperature is assumed). Thus the submillimeter/millimeter 
mass distributions in Perseus and Ophiuchus are quite similar despite different environments, 
distances (d = 160 pc for Ophiuchus, d = 250 pc for Perseus), and resolutions (14" vs. 31"). 
Comparisons are complicated, however, by blending, differing analysis and assumptions, and 
dust property uncertainties. 

We also fit a lognormal distribution to the M > 0.8 M region where the mass function 
is reasonably complete: 
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dM (Znl0)M dlogM' 
Here A is a normalization factor, a is the width of the distribution, and M is the characteris- 
tic mass. The best-fitting lognormal for M > 0.8 M has a = 0.5±0.1 and M = 0.9±0.4 M . 
A lognormal is a somewhat better fit than a broken power law (x 2 = 0.2). For Td = 20 K, 
the best fit parameters are a = 0.5 and M = 0.3 M (note again that only the characteristic 
mass, and not the shape, changes with Td). 

Uncertainties in the dust temperature Td should have only a linear effect on the esti- 
mated mass as long as we are in the Rayleigh- Jeans (RJ) regime of the SED, where B v scales 
linearly with Td- At low temperatures, however, the departure of RJ from a true blackbody 
can cause large errors in the estimated mass. If, for example, the true dust temperature of 
a source is T D = 10 K, assuming a temperature of T D = 5 K would cause a miscalculation 
of the mass by a factor of 4.7, and assuming a temperature of T D = 20 K would result in a 
factor of 3 error. 

Figure 17 demonstrates the effect on the differential mass function of varying the dust 
temperature. For each curve the mass function is calculated for a single value for T D . Note 
that the shape of dN/dM does not change with the assumed value of T D , the distribution 
simply shifts to higher masses (lower T D ) or lower masses (higher T D ). The shape will 
change, however, if there is a range of dust temperatures in the real cores. Ideally we 
could determine T D independently using observations at two different wavelengths. For our 
purposes Td = 10 K is a good compromise for prestellar and protostellar sources, but will 
overestimate by a factor of 3 the mass of sources with a true dust temperature of 20 K. 



Variations in K Llmm and d have similar effects. The dust opacity is uncertain by up 
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to a factor of two or more (Ossenkopf & Henning 1994), owing to large uncertainties in 
the assumed dust properties as well as the possibility that k v varies with position within a 
core. Increasing the opacity shifts the mass distribution to lower masses, while increasing 
the distance shifts it to higher masses. Both are smaller effects than changing the dust 
temperature for the range of plausible values (fti.imm = 0.005—0.02 cm 2 g -1 ; d = 200—300 pc; 
Td — 5 — 30 K). As for variable dust properties within the sample, blending of close sources 
will also distort the shape of the mass function, biasing it toward higher masses. We know 
from previous observations that some sources are blends, but do not attempt to distinguish 
blends from single sources here. The total uncertainties in all masses are at least a factor of 
4 or more (e.g. Shirley et al. 2002; Young et al. 2003). 

A comparison of the prestellar clump mass function to the stellar initial mass function 
(IMF) may reveal the origin of the IMF shape. If stellar masses are determined by compet- 
itive accretion or by the protostars themselves through feedback mechanisms (e.g. outflows 
and winds), we would not expect the emergent IMF to reflect the original clump mass func- 
tion (Adams & Fatuzzo 1996). If, on the other hand, stellar masses are determined by the 
initial fragmentation into cores, as might be expected in crowded regions where the mass 
reservoir is limited to a protostar's nascent core, the IMF should closely trace the clump MF 
(Myers 1998). In Serpens (Testi & Sargent 1998) and p Oph (Motte et al. 1998) the clump 
MF has been found to be quite similar to the stellar IMF, suggesting that fragmentation is 
responsible for determining final masses. In reality, of course, it may well be a combination 
of these effects and turbulence that shapes the IMF (e.g. Clark & Bonnell 2005), but we 
consider the simplest cases here. 

The local IMF follows a broken power law with cti = 1.6 (M < 1 M ) and a 2 = 2.7 
(M > 1 M ), flattening around 0.3 M Q (Chabrier 2003). Chabrier (2003) also find that 
the IMF is well fit by a lognormal with a = 0.4 and M = 0.5 M . Thus the slope of 
the IMF, but not the break mass, is very similar to the Perseus mass function. Even the 
characteristic masses are quite similar if we assume To = 20 K for the Perseus sample. 
Currently, a direct connection between the mass distribution in Perseus and the IMF is 
difficult to make because our sample contains sources at a range of ages, with varying amounts 
of the envelope already accreted onto the protostar or ejected in an outflow, and varying 
envelope temperatures. After combining these data with Spitzer c2d data it will be possible 
to determine the evolutionary state of each source, separating prestellar cores from more 
evolved objects. Even for a sample containing only prestellar cores, however, the association 
with final stellar masses may be problematic (see Johnstone et al. 2000). 
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5.3. Clustering 

Visually, the 1.1 mm sources in Perseus appear very clustered: 89/122 or 73%, have a 
neighboring source within 3', with most isolated sources being faint objects near the detection 
limit. In their 850 /iin SCUBA survey Hatchell et al. (2005) find that 80% of the 850 jum 
sources are in groups of three or more sources (within 0.5 pc), and 40 — 60% of the sources 
are in the massive clusters of NGC 1333 and IC 348 (HH 211 in their paper). Using similar 
criteria, we also find ~ 80% of the total number of sources lie in groups of > 3 within 0.5 pc, 
and ~ 45 — 50% in the massive clusters. Despite the fact that our 1.1 mm map covers more 
than twice as much area as the 850 /im map of Hatchell et al. (2005) (7.5 deg 2 compared 
to ~ 3 deg 2 ), the clustering properties of sources in the two surveys are quite similar. This 
result is perhaps not surprising considering that we only detect about 5 — 10 sources in the 
additional ~ 4.5 deg 2 covered by the Bolocam map. 

For a more quantitative understanding of the clustering properties in Perseus, we cal- 
culate the two-point correlation function: 

where H s (r) is the number of core pairs with separation between log(r) and log(r) +dlog(r). 
The definition of H r (r) is similar to H s (r), but for a random distribution. The random 
sample is constructed by generating a uniform random distribution of sources with the same 
RA and Dec limits as the real sample. The two-point correlation function is often used in 
cosmological studies of the clustering of galaxies (e.g Maddox et al. 1990), but may also be 
a good way to compare the properties of different molecular clouds. 

Plots of H s (r) and w(r) for the entire observed region of Perseus are shown in Figure 18, 
with the random distribution (H r (r), dashed line) included for comparison. If the sources 
were randomly distributed in the cloud, we would expect the two curves to be similar. 
The resolution limit (31" = 8 x 10 3 AU), and average deconvolved source FWHM size 
(1.5 x 10 4 AU) are also indicated, the latter a being representation of the effective resolution 
limit. It is clear that the source pair function H s (r) shows an excess over the random 
distribution H r (r) at small scales (the differences at large scales is not significant). This 
is confirmed by the correlation function w(r) (middle panel), which is > 3a on scales 2 x 
10 4 AU < r < 2.5 x 10 5 AU. Note that the random distribution shows no correlation (w r ~ 0), 
as expected. 

If we characterize the correlation function as a power law, w(r) oc r -7 , then a good fit 
(X 2 = 0.7) is obtained for 7 = 1.25 ± 0.06 in the range 2 x 10 4 AU < r < 2 x 10 5 AU. Since 
the average deconvolved source size corresponds to ~ 1.5 x 10 4 AU, Perseus essentially shows 
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clustering from the average source size up to r max = 2 x 10 5 AU = 1.2 pc. In Ophiuchus 
Johnstone et al. (2000) find 7 = 0.75 for separations r < 3 x 10 4 AU and negligible clustering 
for r > 3 x 10 4 AU, for a distance of 160 pc and a beam size of 2 x 10 3 AU. Those authors 
associate r max = 3 x 10 4 AU with the Jeans length in Ophiuchus, and note that 7 = 0.7 for 
galaxy clustering (Maddox et al. 1990), where gravity is the dominant process. It is possible 
that different slopes could be associated with different processes dominating fragmentation, 
but given the considerable uncertainties involved we choose not to speculate further. 

While the correlation function provides a good comparison of clustering between the 
clouds, we note that the range over which w(r) can be computed (8 x 10 3 — 5 x 10 6 AU for 
Perseus and 4 x 10 3 — 10 5 AU for Oph) differs due to the different resolution limits, spatial 
coverages, and distances of the two clouds. Even so, it seems that clustering occurs on larger 
scales in Perseus than in Oph (r max ,p e r > r maXt o P h), but drops off more quickly as a function 
of separation (7p er > Joph)- It is important to note that the measured correlation function 
depends on the area observed; if we break the Perseus map into smaller pieces the derived 
slope does vary, but it remains between 1 < 7 < 1.5 for a range of chosen areas. 



5.4. An extinction threshold for 1.1 mm cores 

The total mass contained in the 122 detected 1.1 mm cores is 285 M , assuming K, 1Amm = 
0.0114 cm 2 g -1 , T D = 10 K, and a gas to dust ratio of 100. The total cloud mass based on 
CO observations is 1 — 2 x 10 4 M for d = 250 pc (Sargent 1979; Cernicharo et al. 1985; 
Ungerechts & Thaddeus 1987; Carpenter 2000). Based on the NICE extinction map, we 
calculated a total mass for Ay > 2 of 59OOM in the area observed by Bolocam. Thus 
only a very small fraction, between 1% and 5%, of the cloud mass is contained in dense 
cores at A = 1.1mm. This small fraction is consistent with the evidence from comparisons 
of molecular cloud masses to total stellar masses that molecular cloud material is relatively 
sterile (Evans 1999). Hatchell et al. (2005) find a significantly larger fraction of mass (20%) 
in 850 /iin emission. Their large total mass at 850 /im of ~ 2600 M is due in part to 
differences in the assumed distance and opacity, and in part to the fact that those authors 
integrate over all > 5<x pixels in the map to get the total mass, whereas we only include only 
the mass in discrete cores. 

The low efficiency of 1.1 mm cores may be related to the conditions required for a dense 
core to form. Johnstone et al. (2004) have recently suggested that there is an extinction 
threshold for forming 850 /im cores in Ophiuchus. The derived threshold for the presence 
of stable cores in Ophiuchus (A v p l ! ^ m ~ 15 mag) is greater than the maximum extinction 
in Perseus as traced by our NICE extinction map {Ay^ nax ~ 16 mag), suggesting that 
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if conditions in the two regions were similar there should be very few millimeter cores in 
Perseus, which is clearly not the case. In fact, the mean extinction from our NICE map 
toward all cores in Perseus is (Ay) = 7.1 ±3.1 mag, similar to the minimum Ay at which any 
850 /xm core is found in Ophiuchus {Ay P ^ in = 7 mag; Johnstone et al. 2004). Bear in mind, 
however, if the cloud material is very clumpy, beam effects will be important given that the 
distance to Perseus is twice that of Ophiuchus. For comparison, the mean Ay for the entire 
area observed by Bolocam is (A v ) ~ 2 mag, and the mean Ay for the area observed by 
Johnstone et al. (2004) in Ophiuchus was (Ay) ~ 4 mag. 

All dense groups of 1.1 mm sources in Perseus lie in regions of high Ay > 5 mag, 
but not all regions of high A v correspond to 1.1 mm sources, suggesting that relatively high 
extinction may be necessary, but not sufficient, for star formation. The probability of finding 
a 1.1 mm core as a function of Ay is shown in Figure 19. The probability for a given Ay 
is calculated from the extinction map as the number of 50" pixels containing a 1.1 mm core 
divided by the total number of pixels at that Ay (jo/100 = N stc /Na v ). Error bars are Poisson 
statistical errors (cr p /100 = \/N src /NA v )- It appears that there is an approximate extinction 
limit of Ay ~ 5 mag below which it becomes very unlikely that a 1.1 mm core will be found. 
Above Ay ~ 5 mag, the probability of finding a 1.1 mm core rises with Ay. Error bars are 
large for high A v because our NICE Ay map is not sensitive to Ay > 10 mag, so there are 
few pixels at high extinctions. 

Hatchell et al. (2005) take a more sophisticated approach to the probability of finding 
850 fim cores as a function of column density, concluding that there is no column density (or 
Ay) limit for submillimeter cores in Perseus. This result is consistent with our extinction 
limit of A v ~ 5 mag given that the Hatchell et al. (2005) survey covers only the Ay > 4 mag 
regions in Perseus. An extinction threshold around Ay ~ 5 mag in Perseus would also 
explain why very few (5 — 10) sources are found in the Bolocam 1.1 mm map outside of the 
region covered by the Hatchell et al. (2005) survey, despite the much greater area imaged. 
The additional area covered is primarily low column density Ay ~ 2 — 4 mag material, and 
therefore would not be expected to contain many millimeter cores. 

5.5. Comparison to c2d Observations: Bl Ridge 

As an example of the analysis of the complementary c2d Spitzer IR and Bolocam 
1.1 mm observations, we compare the Spitzer 24 /iin and 1.1 mm images for a small area 
of Perseus around Bl. The Bl Ridge is a narrow ridge of extended 1.1 mm emission below 
the group of protostars Bla-d. Figure 20 shows the Bl Ridge region of the Bolocam map 
(left), with the position of all detected MIPS 24 /xm sources, as well as MIPS sources with 
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^24 > 5 mJy, indicated. The c2d MIPS 24 jj.ro. image is also shown (right) with Bolocam 
1.1 mm contours overlaid. The Bl(a-d) protostellar sources are bright at both 24 /im and 
1.1 mm, with (S 24 , S lAmm ) = (0.2,1.2) Jy/beam (Bl-b), (0.8,1.1) Jy/beam (Bl-d) and 
(0.13,0.6) Jy/beam (Bl-d). The little-known protostar IRAS 03292+3039 is relatively faint 
at 24 /im (0.08 Jy/beam), but very bright at 1.1 mm (1.1 Jy/beam), suggesting a young 
evolutionary state. 

There are a few bright 24/xm protostellar sources nearby and in the lower part of the 
extended Bl ridge (e.g. near the ammonia core Per 7 (Ladd et al. 1994)), but in the main 
part of the ridge there are no 24 /im sources. Given the lack of mid-IR sources in the main 
Bl ridge, we suggest that it is made up of a number of prestellar cores. High resolution 
submillimeter and millimeter follow-up observations of both the ridge and nearby protostars 
have already been completed with OVRO (A = 3mm) and SHARC II (A = 350 /im). A 
detailed paper focusing on this region that will combine IRAC and MIPS data with high 
resolution millimeter data to study the precise nature of the prestellar cores and embedded 
protostars, including envelope structure and outflow dynamics, is in preparation. 



6. Summary 

We present a 7.5 deg 2 (143 pc 2 for d = 250 pc) survey for 1.1 mm dust continuum 
emission in the Perseus molecular cloud using Bolocam at the CSO. This map is the largest 
millimeter or submillimeter continuum image of the region to date. Given that Bolocam has 
a beam of 31" (FWHM), the Perseus map covers a remarkable total number (3.4 x 10 4 ) of 
resolution elements down to a 5a point source mass detection limit of 0.18 M Q (T D = 10 K). 

We detect 122 1.1 mm cores above the 5a = 75mJy/beam detection limit. Nearly half 
(60/122) of the detected sources are new millimeter or submillimeter detections, including 
previously unknown sources as well as known objects not previously observed at these wave- 
lengths. Our 1.1 mm map covers more than twice the area of the recent 850 /xm SCUBA 
survey by Hatchell et al. (2005) but, despite the significantly greater area imaged, only 
~ 5 — 10 of our new detections lie outside their map. Thus, much of the Perseus cloud is 
devoid of compact millimeter emission and, by implication, active star formation. 

In order to compare the general morphology of our 1.1 mm map to the COMPLETE 
13 CO map and our NICE extinction map, we convert all three images to a column density 
scale. Our 1.1 mm map reveals significantly higher column density features than the other 
tracers and exhibits much more compact structure, even when degraded to the 5' resolution 
of the NICE extinction map. The general appearance of the 1.1 mm emission is roughly 
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consistent with the molecular and extinction data, however, in that most 1.1 mm sources lie 
within 13 CO and Ay peaks. 

The total mass in discrete 1.1 mm cores is 285 M Q (T D = 10 K), accounting for no more 
than 5% of the total mass of the cloud. The small fraction of mass in dense cores, which are 
usually associated with star formation, supports the idea that most of the mass in molecular 
clouds is relatively sterile (Evans 1999). Calculating the probability of finding a 1.1 mm core 
as a function of Ay leads us to conclude that there is an extinction threshold in Perseus at 
Ay ~ 5 mag, above which 1.1 mm cores are likely to be observed. Such an extinction limit 
is consistent with the fact that very few new sources are found outside of the area covered 
by the Ay > 4 mag map of Hatchell et al. (2005). 

The average mass of the sample, based on the total flux density from aperture photome- 
try, is 2.3 M & (Tq = 10 K). The differential mass function dN/dM is well fitted by a broken 
power law with a ~ 1.3 (0.5 M. < M < 2.5M ) and a ~ 2.6 (M > 2.5 M ). The derived 
values are similar to those found in Ophiuchus (a = 2.5 M > 1 M Q ,a — 1.5 M < 1 M Q 
Johnstone et al. 2000), and to the local IMF (a = 2.7 M > 1 M Q ,a = 1.6 M < 1 M , 
Chabrier 2003). We compute the two-point correlation function, confirming that compact 
millimeter emission in the cloud is highly structured. Significant clustering of sources from 
the average source size up to scales of 2 x 10 5 AU is seen. Within this range, the correlation 
function exhibits a power law shape with index 7 = —1.25, steeper than the correlation in 
Ophiuchus (Johnstone et al. 2000). 

The Bolocam 1.1 mm data presented here were designed to cover the same region as 
the Spitzer c2d legacy program IRAC (A = 3.6 - 8.0 /im) and MIPS (A = 24 - 160 /iin) 
maps of Perseus. Combining the 1.1 mm data with the c2d IRAC and MIPS data will en- 
able a complete census of the properties and distribution of the protostars and dense cores 
in Perseus, especially when complemented by the publicly available 2MASS catalogs and 
COMPLETE molecular and continuum maps. Follow-up observations have already been 
completed with SHARC II (CSO, A = fim) and the OVRO interferometer (A = 3mm) for 
several of the most interesting sources identified in the Bolocam map, including the Bl ridge. 
A detailed discussion of these regions utilizing Spitzer photometry, high resolution submil- 
limeter and millimeter images, molecular line data, and SED and radial profile modeling is 
in preparation. 
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Fig. 1. — Bolocam 1.1 mm (thin line) and Spitzer c2d IRAC (thick line) coverage of Perseus 
overlaid on a 13 CO integrated intensity map from Padoan et al. (1999). The area observed 
by IRAC was chosen based on the 13 CO intensity, and corresponds approximately to Ay > 
2 mag (Evans et al. 2003). The Bolocam observations were designed to cover the same 
region. The MIPS data cover a somewhat larger area. 
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Fig. 2. — NGC 1333, before (left) and after (center) iterative mapping, shown on the same 
intensity scale. Dark blue regions are negative lobes introduced by cleaning. These im- 
ages illustrate the effectiveness of iterative mapping in reducing such negative artifacts and 
restoring source flux density lost during PCA cleaning. After 10 iterations, the brightest 
source has increased in peak brightness by 14%. Although some negative artifacts remain, 
they are greatly reduced both in extent and intensity (the most negative pixel has decreased 
in amplitude from —238 mJy/beam to —88 mJy/beam). No fine-scale structure is lost, and 
all recovered extended structure is real; for comparison we show the 850 /jm map of Sandell 
& Knee (2001), with a resolution of 14", on approximately the same intensity scale. 
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input S/N 

Fig. 3. — Fractional lost peak flux density of simulated sources after 10 iterations, as a 
function of the input signal to noise (S/N) and FWHM. Solid lines are the 1 sigma RMS 
noise as a fraction of the input amplitude, indicating the spread in recovered peak that might 
be expected from noise alone. Except for the largest (120" FWHM) sources, most points 
lie within the RMS noise, so residual systematic effects from cleaning are not important for 
measured peak flux densities. 
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Fig. 4. — Fractional lost integrated flux density in a 40" aperture for simulated sources, as 
a function of input signal to noise (S/N), FWHM, and iteration number. Representative 
source FWHM of 30", 50", 80", 120" are shown, with insets to magnify confused regions. 
Larger sources require more iterations to recover the input flux density, and tend to have 
the most remaining fractional missing flux density after 10 iterations. Except for sources 
with FWHM> 100" the flux has converged by ~ 5 iterations. The integrated flux density 
is usually recovered to within 5%, and to within 10% even for large (FWHM> 80"), faint 
(S/N< 20) sources. Note that the lex RMS in a 40" aperture is > 10% for sources with 
S/N< 18. 



-38- 




Fig. 5. — Bolocam 1.1 mm map of the Perseus Molecular Cloud (10"/pixel). The 31" 
resolution map covers 7.5 deg 2 (143 pc 2 at d — 250 pc), or 3.4xl0 4 resolution elements. The 
average la RMS is 15 mJy/beam, varying by 15% across the map due to variable coverage. 
In this and other figures all maps shown are unfiltered. 
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Fig. 6. — Bolocam 1.1 mm map of Perseus, with the 122 1.1 mm sources detected above 
5cr indicated by small circles. The RMS varies very little across the map (15%), so the 
apparent lack of sources over large regions of the image is real. Despite the greater area 
surveyed compared to previous work, few new sources are found; much of the cloud is devoid 
of 1.1 mm emission at this sensitivity. Regions of high source density are magnified for 
clarity. 
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Fig. 7. — Examples of new millimeter detections. Each image is 8.5' x 8.5', and has been 
smoothed to an effective resolution of 35". Contours are (2, 4, 6, 8, 12, 16, 20, 25, 35, 50, 75) cr 
on the grayscale 1.1 mm map. Any well know sources are labeled with their common names. 
Numbers give the "Bolo#" identification from this paper (Table 1); those sources with a "*" 
next to the ID were either not covered by or not detected in the 850 /im SCUBA survey of 
Hatchell et al. (2005). Sources range from compact (B) to extended (L) and crowded (I) to 
isolated (C). 
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(0002) 

Fig. 8. — Visual extinction (A v ) contours calculated from 2MASS data using the NICE 
method (see Section 4.2), overlaid on the grayscale 1.1 mm map. Contours in the following 
plots are Ay = 2, 3, 4, 6, 8 mag with an effective resolution of 5' unless otherwise noted. Most 
1.1 mm sources lie within relatively high extinction (Ay > 5) peaks. 
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Fig. 9. — Ay contours calculated from the COMPLETE 13 CO map, smoothed to 5' to match 
the resolution of the NICE extinction map, overlaid on the Bolocam 1.1 mm map. Thick 
black lines indicate the observational boundaries of the COMPLETE map. As a tracer of 
column density, 13 CO is only accurate where it is optically thin. A v ( 13 CO) becomes more 
compact compared to extinction Ay, especially around 1.1 mm sources. 
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Fig. 10. — Ay contours calculated from the 1.1 mm map, smoothed to 5' resolution. Contours 
are Ay = 2, 5, 10, 15, 20 mag. The 1.1 mm contours are much more compact and extend to 
higher column density than the other column density tracers, likely because 1.1 mm emission 
traces only the densest regions of the cloud. 
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Fig. 11. — Distribution of source peak flux densities (dashed line, m Jy/beam) and total flux 
densities (solid line, Jy) from aperture photometry. Clumpfind 3a flux densities are also 
shown for comparison (dash-dot line, Jy). The vertical dotted line is the 5cr peak detection 
limit. The peak flux density distribution has a mean of 0.35 Jy/beam and dispersion of 
0.56 Jy/beam (0.35 ± 0.56), and the total flux density distribution a mean of 0.96 ± 1.21 Jy. 
The mean flux from Clumpfind (1.13 ± 1.80 Jy) is somewhat larger than that from aperture 
photometry. 
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Fig. 12. — Distribution of minor axis (dashed line) and major axis (solid line) FWHM sizes, 
as determined from an elliptical gaussian fit. The beam size is indicated by the dotted line. 
The mean minor axis FWHM is 58" ± 17", and the mean major axis FWHM is 80" ± 27". 
Full-width at 3a sizes for Clumpfmd sources are also shown (dash-dot line), where (FW3a = 
2 y/ Npi x Api X 1 7T . The mean size from Clumpfmd (77" ± 30") is larger than the mean FWHM 
(68" ± 20"), as expected because the Clumpfmd size is measured at the 3a contour rather 
than at the half-max. 
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Fig. 13. — Distribution of source axis ratios. The mean axis ratio is 1.4, and the sample 
contains some very elongated sources with axis ratios > 2. We find that measured axis ratios 
< 1.2 are unreliable based on Monte Carlo tests. 
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Fig. 14. — Total mass vs. FWHM size for Peak-find sources. The solid line denotes the 
analytic mass detection limit as a function of size (M oc R 2 ) for gaussian sources. Empir- 
ical 10%, 50%, and 90% completeness curves are also shown, derived using Monte Carlos 
with simulated sources and taking into account the effects of cleaning, iterative mapping, 
and optimal filtering. 50% completeness for Bonnor-Ebert spheres (thick shaded line) were 
similarly determined using BE models with n c = (8 x 10 4 , 8 x 10 4 , 9.5 x 10 4 ) cm" 3 and 
r a = (8 x 10 3 , 1.5 x 10 4 , 3 x 10 4 ) AU. Representative error bars for 50" and 100" FWHM 
sources near the detection limit are shown, as estimated from the results of Monte Carlo 
simulations. Note the lack of sources near the resolution limit, which cannot be entirely be 
accounted for by pointing errors of < 7", although the pointing-smeared beam could be as 
large as 34" (shaded region). 
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Fig. 15. — Total mass vs. FW3cr size for Peak-find and Clumpfind sources. The difference 
in the two distributions lies in the total flux calculation: aperture photometry for Peak-find, 
and integrating over the > 3a pixels for Clumpfind. Sizes for Peak-find sources are computed 
by scaling the FWHM to the 3a contour (FW3cr =FWHMx ^\n(S U:Peak /3a) / ln2), assuming 
the source is gaussian. The Clumpfind size is based on the total pixel area above the 3a 
contour (FW3<r = 2^Npi X A pix /iT). Completeness curves are as in Figure 14, scaled to 
the 3a size. The smaller apparent scatter compared to Figure 14 is not real, but rather a 
consequence of the different size definitions used. 
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Fig. 16. — Differential mass function dN/dM for masses calculated from aperture photom- 
etry and a single dust temperature To = 10 K. The 50% completeness limit is 0.18 M 
{To = 10 K) for a point source, and 0.8 M for a 70" FWHM source, the average size of 
the sample. Assuming a broken power law of the form N(M) oc M~ a , the best fit slopes are 
ai = 1.3 ±0.3 (O.5M < M < 2.5M ) and a 2 = 2.6 ±0.3 (M > 2.5 M Q ). The slope, but not 
the break mass, is very similar to the local IMF: a± = 1.6 (M < 1 M ), a 2 = 2.7 (M > 1 M & ) 
(Chabrier 2003). The data for M > 0.8 M , where the distribution is not affected by com- 
pleteness, is also well fitted by a lognormal with a = 0.5 ±0.1, M = 0.9 ±0.4 M . The best 
fit single power law is a = 2.1 ± 0.1 for M > 0.5 M . For comparison, 850 /zm sources in 
Ophiuchus were found to have a\ ~ 1.5 below 1 M and a 2 ~ 2.5 above 1 M (Johnstone 
et al. 2000). 
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Fig. 17. — The differential mass function, shown for different values of the assumed dust 
temperature Tp. The dashed lines correspond to the mass limit for a point source at each T/> 
Changing the temperature shifts the distribution to higher (for lower Tp) or lower (for higher 
Td) masses, but does not change its shape. A dust opacity of «i.i mm = 0.0114 cm 2 g _1 and a 
distance of d — 250 pc are assumed for all masses. Increasing K\.x mm shifts the distribution 
to lower masses, while increasing d shifts it to higher masses. Given the range of plausible 
values («4.i mm = 0.005 — 0.02 cm 2 /g, d = 200 — 300 pc), the effects of varying «4.i mm and d 
are smaller than the effect of varying the temperature (Td = 5 — 30 K). 
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Fig. 18. — The two-point correlation function, illustrating the degree of clustering in the 
cloud. Top: H(r) is the number of source pairs between log(r) and dlog(r) as a function 
of log(r/AU). The solid line indicates the real data, and the dashed line is for a uniform 
random distribution of sources with the same RA/Dec limits as the real sample. In all plots 
the dotted line denotes the resolution limit, and the dot-dash line the mean deconvolved 
source FWHM. Center: The two point correlation function w(r) as defined in the text, with 
VN errors. Where w(r) > there is a correlation between sources at that separation, thus 
w(r) indicates clustering on scales log(r / AU) < 5.5 for the Perseus sample. When calculated 
using two randomly distributed samples (dashed line) w(r) shows no correlation, as expected. 
Bottom: w(r) is well fitted by a power law, w(r) oc r _L25 , for 2.5 x 10 4 AU < r < 2 x 10 5 AU. 
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Fig. 19. — Probability of finding a 1.1 mm core as a function of A v , computed using the 
NICE extinction map. The probability is the number of 50" pixels at a given Ay containing 
one or more 1.1 mm cores, divided by the total number of pixels at that Ay. Error bars are 
Poisson statistical errors, and become large where the Ay map saturates at Ay ~ 10 mag. 
The low Ay region is magnified for clarity. It appears that there is an approximate extinction 
limit at Ay ~ 5 mag, below which it becomes very unlikely that a 1.1 mm core will be found. 
The extinction limit of Ay ~ 5 mag is considerably lower than the Ay ~ 15 mag limit for 
forming 850 /im cores found by Johnstone et al. (2004) in Ophiuchus, although beam effects 
could be important. Our Ay ~ 5 mag limit is consistent with the fact that few sources are 
found outside previously surveyed regions. 
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Fig. 20. — Spitzer sources in the Bl Ridge region overlaid on the Bolocam 1.1 mm map (left). 
Blue diamonds show the positions of all detected 24 /j,m sources in the MIPS c2d image 
(Rebull et al. 2005, in preparation), and pink boxes indicate sources with S24 > 5 mJy. The 
MIPS 24 /iin image is also shown (right) with Bolocam 1.1 mm contours overlaid (2, 4, ..20er). 
There are no 24 /zm sources in the main part of the Bl Ridge, suggesting that the ridge may 
be made up of a number of prestellar cores. 
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Table 1. Source IDs, positions, and peak flux density 



ID 


RA (2000) 
(h m s) 


Dec (2000) 
(° ' ") 


Peak 
(mjy/bcam) 


S/N 


Other Names 


MIPS 24^m Source? 


Bolo 1 


03 25 07.8 


+30 24 21.6 


109 


(16) 


6.1 




N 


Bolo 2 


03 25 09.5 


+30 23 51.3 


121 


(16) 


8.7 




N 


Bolo 3 


03 25 10.2 


+30 44 43.4 


125 


(17) 


6.6 




N 


Bolo 4 


03 25 17.0 


+30 18 53.2 


149 


(17) 


7.2 




N 


Bolo 5 


03 25 22.2 


+30 45 09.3 


727 


(18) 


42.0 


L 1448-IRS2 (1) 


Y 


Bolo 6 


03 25 26.5 


+30 21 50.0 


143 


(19) 


5.9 




N 


Bolo 7 


03 25 35.4 


+30 13 06.2 


126 


(17) 


6.0 




N 


Bolo 8 


03 25 35.9 


+30 45 17.2 


2241 


(19) 


128 


L 1448-N (2) 


Y 


Bolo 9 


03 25 37.2 


+30 09 55.4 


113 


(16) 


5.7 




N 


Bolo 10 


03 25 38.4 


+30 43 58.0 


980 


(20) 


45.2 


L 1448-C (2) 


Y 


Bolo 11 


03 25 46.0 


+30 44 10.3 


240 


(20) 


8.4 




Y 


Bolo 12 


03 25 47.5 


+30 12 25.9 


118 


(18) 


5.1 




N 


Bolo 13 


03 25 48.8 


+30 42 24.1 


407 


(19) 


19.0 


H32 (3) 


N 


Bolo 14 


03 25 50.6 


+30 42 01.4 


342 


(19) 


17.9 


H33 (3) 


N 


Bolo 15 


03 25 55.1 


+30 41 25.8 


188 


(17) 


6.4 




N 


Bolo 16 


03 25 56.2 


+30 40 41.1 


154 


(17) 


6.3 




N 


Bolo 17 


03 25 58.4 


+30 37 13.4 


122 


(15) 


6.7 




N 


Bolo 18 


03 26 36.9 


+30 15 23.0 


225 


(15) 


14.9 


IRAS 03235+3004; H80 (3) 


Y 


Bolo 19 


03 27 02.0 


+30 15 07.9 


115 


(14) 


7.3 




N 


Bolo 20 


03 27 28.8 


+30 15 02.1 


114 


(16) 


5.9 




N 


Bolo 21 


03 27 37.5 


+30 13 53.4 


164 


(17) 


9.0 


H39 (3) 


Y 


Bolo 22 


03 27 39.3 


+30 12 52.9 


317 


(18) 


16.3 


RNO 15-FIR (4) 


Y 


Bolo 23 


03 27 41.7 


+30 12 24.1 


283 


(18) 


12.8 


H36 (3) 


Y 


Bolo 24 


03 27 47.9 


+30 12 02.1 


220 


(19) 


11.2 


L 1455-IRS2 (5); H37 (3) 


Y 


Bolo 25 


03 28 32.1 


+31 11 08.8 


106 


(16) 


6.7 


H74 (3) 


N/A 


Bolo 26 


03 28 32.4 


+31 04 42.9 


117 


(14) 


7.6 




N/A 


Bolo 27 


03 28 33.3 


+30 19 35.0 


155 


(16) 


6.4 




Y 


Bolo 28 


03 28 34.0 


+31 07 01.2 


137 


(14) 


9.7 


H69 (3) 


N/A 


Bolo 29 


03 28 36.2 


+31 13 26.4 


267 


(19) 


17.1 


NGC 1333-IRAS 1 (6); SK6 (7) 


N/A (IRAS) 


Bolo 30 


03 28 39.0 


+31 05 59.2 


189 


(15) 


11.3 


in Per3 (8); H71 (4) 


N/A 


Bolo 31 


03 28 39.9 


+31 17 56.7 


385 


(19) 


20.1 


NGC 1333-IRAS 5 (6) 


N/A (IRAS) 


Bolo 32 


03 28 41.7 


+30 31 13.0 


144 


(15) 


6.4 


RNO 17 


Y 


Bolo 33 


03 28 42.6 


+31 06 12.2 


173 


(14) 


9.6 


in Per 3 (8); H75 (3) 


N/A 


Bolo 34 


03 28 45.9 


+31 15 19.8 


150 


(20) 


6.2 




N/A 


Bolo 35 


03 28 48.4 


+31 16 01.9 


130 


(20) 


5.9 




N/A 


Bolo 36 


03 28 48.7 


+30 43 24.6 


127 


(14) 


6.3 




Y 


Bolo 37 


03 28 52.1 


+31 18 07.9 


160 


(20) 


7.4 


ASR 40 


N/A 


Bolo 38 


03 28 55.3 


+31 14 32.2 


1220 


(20) 


59.4 


NGC 1333-IRAS 2 (6,9) 


N/A (IRAS) 


Bolo 39 


03 28 55.3 


+31 19 17.7 


210 


(20) 


9.7 


ASR 64; H68 (3) 


N/A 


Bolo 40 


03 29 00.0 


+31 21 37.8 


685 


(17) 


34.2 


SK31 (7) 


N/A 


Bolo 41 


03 29 00.6 


+31 11 58.9 


190 


(19) 


10.7 


SKI (7) 


N/A 


Bolo 42 


03 29 01.3 


+31 20 33.0 


1079 


(18) 


58.8 


NGC 1333-IRAS 6 (6); SK24 (7) 


N/A (IRAS) 


Bolo 43 


03 29 02.6 


+31 15 56.8 


2160 


(20) 


103 


SVS 13 (9) 


N/A 


Bolo 44 


03 29 04.5 


+31 18 42.1 


270 


(20) 


13.7 




N/A 


Bolo 45 


03 29 07.7 


+31 17 16.8 


460 


(20) 


21.8 


Cor 1 (9); SK17/18 (7) 


N/A 
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Table 1 — Continued 



ID RA (2000) Dec (2000) Peak S/N 







(h 


m s 


) 


(° 
\ 




"1 




(mjy/beam) 




Bolo 


46 


03 


on 
Z\) 


OQ 


3 


+31 


15 


11. 


1 


740 (20) 


28.6 


Bolo 


47 


03 


on 
Z\) 


Uo 


,9 


+31 


21 


43. 


5 


612 (18) 


30.9 


Bolo 


48 


03 


oo 
Z\) 


1 o 
1U. 


.9 


+31 


13 


26 


4 


5180 


(20) 


264 


Bolo 


49 


03 


Z\) 




3 


+31 


18 


24. 


.8 


840 


(20) 


40.6 


Bolo 


50 


03 


OO 
Z\) 




,5 


+31 


20 


30. 


1 


310 


(20) 


14.2 


Bolo 


51 


03 


Z\) 


1 7 
1 1 





+31 


12 


25. 


9 


423 


(19) 


16.4 


Bolo 


52 


03 


on 
Z\) 


1 / 


,1 


+31 


27 


38. 


8 


269 


(18) 


13.7 


Bolo 


53 


03 


oo 

Z\) 


ICS. 


.2 


+31 


25 


17. 





336 


(19) 


16.7 


Bolo 


54 


03 


oo 
z\> 


1 o 

ly. 


,1 


+31 


23 


26. 


5 


330 


(19) 


16.7 


Bolo 


55 


03 


on 
Z\) 


111 


,4 


+31 


11 


36 


6 


184 


(17) 


9.8 


Bolo 


56 


03 


Z\) 


oo 
zz 


,4 


+31 


36 


24. 


.4 


102 


(14) 


6.8 


Bolo 


57 


03 


Z\> 


00 

zz 


.9 


+31 


33 


16 


5 


224 


(16) 


16.2 


Bolo 


58 


03 


90 
Z\) 


zo. 


.7 


+31 


28 


16. 


3 


273 


(17) 


14.0 


Bolo 


59 


03 


90 
Z\) 


ol 


,5 


+31 


39 


12 


9 


249 


(14) 


19.3 


Bolo 


60 


03 


on 
OU 


1 -1 
11 


.9 


+30 


23 


36. 


9 


125 


(15) 


8.3 


Bolo 


61 


03 


on 

OU 


OA 
Z 1 





+30 


27 


38. 


5 


106 


(17) 


5.1 


Bolo 


62 


03 


OU 


QO 
OZ 





+30 


26 


18. 


6 


199 


(17) 


11.1 


Bolo 


63 


03 


on 

OU 


/1 ~, 

40 


,5 


+30 


52 


34. 


3 


141 


(14) 


10.2 


Bolo 


64 


03 


on 
ou 


^n 
ou. 


.4 


+30 


49 


17. 


4 


86 (14) 


6.4 


Bolo 


65 


03 


Q 1 
Ol 


on 
zU. 


.0 


+30 


45 


30. 


2 


522 


(17) 


31.7 


Bolo 


66 


03 


oo 

oz 


1 7 
1 i 


3 


+30 


49 


44 





1068 


(15) 


73.5 


Bolo 


67 


03 


on 

oz 


Oft 
ZD. 


,9 


+30 


59 


11. 





165 


(19) 


8.8 


Bolo 


68 


03 


QO 
OZ 


no 
ZO 


,1 


+31 


02 


17. 


5 


110 


(18) 


6.6 


Bolo 


69 


03 


to 
oz 


QO 

oy 


3 


+30 


57 


28. 


.8 


155 


(16) 


5.7 


Bolo 


70 


03 


OZ 


'ii 


.0 


+31 


00 


00. 





238 


(17) 


11.1 


Bolo 


71 


03 


QO 
OZ 


Ol. 


,2 


+31 


01 


47. 


6 


180 


(20) 


7.3 


Bolo 


72 


03 


32 


57 


.0 


+31 


03 


20 


.8 


240 


(20) 


9.4 


Bolo 


73 


03 


33 


00 





+31 


20 


43. 


8 


178 


(16) 


8.8 


Bolo 


74 


03 


33 


01 


.9 


+31 


04 


31. 


8 


260 


(20) 


9.3 


Bolo 


75 


03 


33 


04 


,4 


+31 


04 


58, 


.8 


260 


(20) 


9.7 


Bolo 


76 


03 


33 


11 


,4 


+31 


21 


30. 


9 


116 


(16) 


6.5 


Bolo 


77 


03 


33 


11 


.4 


+31 


17 


24. 





110 


(19) 


5.0 


Bolo 


78 


03 


33 


13. 


.2 


+31 


19 


50. 


1 


262 


(16) 


16.3 


Bolo 


79 


03 


33 


15. 


.0 


+31 


07 


01. 


9 


605 


(19) 


30.0 


Bolo 


80 


03 


33 


17 


,8 


+31 


09 


29 


8 


1080 


(20) 


53.2 


Bolo 


81 


03 


33 


20 


.5 


+31 


07 


37. 


2 


1220 


(20) 


59.9 


Bolo 


82 


03 


33 


25 


,1 


+31 


05 


34. 


.8 


147 


(17) 


7.1 


Bolo 


83 


03 


33 


25 


.4 


+31 


20 


07. 


4 


165 


(17) 


7.2 


Bolo 


84 


03 


33 


26 


.8 


+31 


06 


49. 


6 


190 


(20) 


7.0 


Bolo 


85 


03 


33 


32 


,1 


+31 


20 


03. 


.8 


155 


(17) 


6.9 


Bolo 


86 


03 


33 


51 


,1 


+31 


12 


36. 


.7 


131 


(18) 


5.2 


Bolo 


87 


03 


35 


21 


.6 


+31 


06 


55. 


4 


99 (14) 


5.9 


Bolo 


88 


03 


40 


14 


,5 


+32 


01 


29. 


6 


162 (17) 


6.7 


Bolo 


89 


03 


40 


49 


3 


+31 


48 


34. 


5 


151 (16) 


8.1 


Bolo 


90 


03 


41 


09 





+31 


44 


33. 





139 (15) 


8.0 



Other Names MIPS 24^m Source? 



SK16 (7) 


N/A 


SK28 (7) 


N/A (IRAS) 


NGC 1333-IRAS 4 (6,10) 


N/A (IRAS) 


NGC 1333-IRAS 7 (6,9) 


N/A (IRAS) 


SK22 (7) 


N/A 


H59 (3) 


N/A 


in Per 4 (8); H61 (3) 


N/A 


SK33 (7) 


N/A (IRAS) 


SK32 (7) 


N/A 


ASR 129; H72 (3) 


N/A 


H91 (3) 


N/A 


in Per 4 (8); H58 (3) 


N/A (IRAS) 


H64 (3) 


N/A 


Per 5 (8) 


N/A (IRAS) 


IRAS 03271 (11); H15 (3); Per 6 (8) 


Y 




N 




Y 




N/A (IRAS) 




N/A 


IRAS 03282 (11); H77 (3) 


N/A (IRAS) 


IRAS 03292+3039; H76 (3) 


Y 


H89 (3) 


N 




Y 




N 


in Per 7 (8) 


N 




N 




N 


in Per 9 (8) 


N 


H5 (3) 


N 


H8 (3) 


N 


in Per 9 (8) 


N 




N 


in Per 9 (8); H82 (3) 


Y 


Bl-d (12) 


Y 


Bl-c (12) 


Y 


Bl-b (12) 


Y 




N 


in Per 9 (8) 


N 


H10 (3) 


Y 


in Per 9 (8) 


N 




Y 




N 




N 




Y 




Y 
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Table 1 — Continued 



ID 


RA (2000) 
(h m s) 


Dec (2000) 
(° ' ") 


Peak 
(mjy/bcam) 


S/N 


Other Names 


MIPS 24^m Source? 


Bolo 91 


03 41 19.7 


+31 47 28.3 


106 (15) 


6.5 




N 


Bolo 92 


03 41 40.0 


+31 58 04.8 


94 (14) 


6.4 




N 


Bolo 93 


03 41 45.1 


+31 48 09.7 


118 (16) 


6.3 




N 


Bolo 94 


03 41 45.8 


+31 57 21.6 


132 (14) 


7.3 




N 


Bolo 95 


03 42 20.6 


+31 44 49.2 


108 (14) 


6.2 




Y 


Bolo 96 


03 42 47.1 


+31 58 40.8 


148 (16) 


7.1 




N 


Bolo 97 


03 42 52.3 


+31 58 12.3 


149 (17) 


6.0 




Y 


Bolo 98 


03 42 57.3 


+31 57 49.3 


134 (18) 


6.0 




Y 


Bolo 99 


03 43 38.3 


+32 03 08.6 


164 (13) 


12.0 


H23 (3) 


N/A 


Bolo 100 


03 43 44.0 


+32 03 10.0 


250 (14) 


15.3 


H26 (3) 


N/A 


Bolo 101 


03 43 45.6 


+32 01 45.1 


139 (19) 


7.5 




N/A 


Bolo 102 


03 43 50.5 


+32 03 17.2 


432 (16) 


23.5 


H81 (3) 


N/A 


Bolo 103 


03 43 55.9 


+32 00 46.4 


994 (19) 


56.5 


HH 211 (13) 


N/A (IRAS) 


Bolo 104 


03 43 57.3 


+32 03 03.9 


777 (17) 


44.2 


IC 348-MMS (14) 


N/A 


Bolo 105 


03 43 57.8 


+32 04 06.2 


283 (16) 


12.5 


H17 (3) 


N/A 


Bolo 106 


03 44 01.6 


+32 02 02.0 


426 (18) 


18.1 


H16 (3) 


N/A 


Bolo 107 


03 44 02.1 


+32 02 33.7 


388 (18) 


16.5 


H18 (3) 


N/A 


Bolo 108 


03 44 02.3 


+32 04 57.3 


139 (16) 


6.1 




N/A 


Bolo 109 


03 44 05.0 


+32 00 27.7 


127 (17) 


6.3 




N/A 


Bolo 110 


03 44 05.2 


+32 02 05.6 


257 (19) 


11.8 


H20 (3) 


N/A 


Bolo 111 


03 44 14.5 


+31 58 00.1 


170 (17) 


8.9 




N/A 


Bolo 112 


03 44 14.8 


+32 09 12.6 


127 (15) 


8.2 




N/A 


Bolo 113 


03 44 22.6 


+31 59 23.2 


122 (18) 


6.2 




N/A 


Bolo 114 


03 44 22.7 


+32 10 01.2 


140 (15) 


8.0 




N/A 


Bolo 115 


03 44 36.4 


+31 58 39.3 


167 (17) 


10.4 


H19 (3) 


N/A 


Bolo 116 


03 44 43.9 


+32 01 24.6 


289 (18) 


18.8 


H14 (3) 


N/A (IRAS) 


Bolo 117 


03 44 48.8 


+32 00 29.5 


150 (18) 


8.7 


H25 (3) 


N/A 


Bolo 118 


03 44 56.0 


+32 00 31.3 


109 (17) 


5.7 




N/A 


Bolo 119 


03 45 15.9 


+32 04 48.3 


181 (16) 


9.3 


H90 (3) 


N/A 


Bolo 120 


03 45 48.1 


+32 24 15.1 


83 (10) 


8.8 


LkHA 330 


Y 


Bolo 121 


03 47 33.5 


+32 50 54.9 


136 (15) 


6.3 




N 


Bolo 122 


03 47 40.8 


+32 51 57.2 


262 (17) 


14.3 


B5-IRS1 (15) 


Y 



Note. — Numbers in parentheses are Icr uncertainties, and do not include an additional systematic un- 
certainty of 15% in the peak flux density. Other names listed are the most common identifications from 
the literature, and are not meant to be a complete list. Presence of MIPS source is from the c2d database 
(Stapelfeldt et al. 2005, in preparation). References - (1) O'Linger et al. 1999; (2) Barsony et al. 1998; (3) 
Hatchell et al. 2005; (4) Froebrich et al. 2003; (5) Juan et al. 1993; (6) Jennings et al. 1987; (7) Sandell & 
Knee 2001; (8) Ladd et al. 1994; (9) Lefloch et al. 1998; (10) Sandell et al. 1991; (11) Bachiller et al. 1991; 
(12) Mathews & Wilson 2002; (13) McGaughrean et al. 1994; (14) Eisloffel et al. 2003; (15) Motte & Andre 
2001. 



Table 2. Source photometry, masses, column density, sizes, and morphology 



ID 


Flux(40") 


Flux(80") 


Flux(120") 


Total Flux 


Mass (10K) Peak A v 


FWHM FWHM 


PA 


<n> 


Morphology 1 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


(M©) (mag) 


(minor,") (major,") 


(°) 


(cm 3 ) 





Bolo 


1 


















0.114 


(0.016) 


0.27 


(0.04) 


8 


38 


(1 


1) 


57 (1.7) 


-52 (6) 


lxlO 5 


multiple, round, weak 


Bolo 


2 


















0.148 


(0.016) 


0.36 


(0.04) 


9 


33 


(1 


o) 


54 (1.6) 


48 (6) 


3xl0 5 


multiple, round, weak 


Bolo 


3 





15 


(0.02) 





30 


(0.05) 


0.38 


(0.07) 


0.38 


(0.07) 


0.9 


(0.2) 


9 


65 


(1 


3) 


134 (3) 


-68 (2) 


3xl0 4 


extended, round, weak 


Bolo 


4 





21 


(0.02) 





48 


(0.05) 


0.87 


(0.07) 


1.24 


(0.09) 


3.0 


(0.2) 


11 


99 


(1 


2) 


122 (1.4) 


74 (4) 


5xl0 4 


extended, round 


Bolo 


5 





93 


(0.02) 


1 


48 


(0.05) 


2.05 


(0.07) 


2.05 


(0.07) 


4.9 


(0.2) 


51 


45 


(0 


2) 


79 (0.4) 


88 (1) 


7xl0 5 


multiple, extended, elongated 


Bolo 


6 





19 


(0.03) 





43 


(0.05) 


0.79 


(0.08) 


1.18 


(0.10) 


2.9 


(0.2) 


10 


93 


(1 


3) 


132 (1.8) 


-54 (3) 


4xl0 4 


extended, round, weak 


Bolo 


7 





16 


(0.02) 





36 


(0.04) 


0.67 


(0.07) 


0.96 


(0.09) 


2.3 


(0.2) 


9 


88 


(1 


4) 


122 (1.9) 


4 (4) 


4xl0 4 


multiple, extended, round, weak 


Bolo 


8 


2 


95 


(0.03) 


4 


70 


(0.05) 






4.70 


(0.05) 


11.3 


(0.12) 


160 


45 


(0 


1) 


54 (0.1) 


-47 (1) 


4xl0 6 


multiple, extended, round 


Bolo 


9 





15 


(0.02) 





31 


(0.04) 


0.52 


(0.07) 


0.70 


(0.09) 


1.7 


(0.2) 


8 


76 


(1 


3) 


144 (3) 


-52 (2) 


3xl0 4 


extended, round, weak 


Bolo 


10 


1 


22 


(0.03) 


2 


01 


(0.06) 






2.04 


(0.06) 


4.9 


(0.14) 


69 


53 


(0 


2) 


55 (0.2) 


-10 (10) 


lxlO 6 


multiple, extended, round 


Bolo 


11 





31 


(0.03) 





75 


(0.05) 






0.90 


(0.06) 


2.2 


(0.15) 


17 


78 


(0 


8) 


81 (0.8) 


30 (20) 


lxlO 5 


multiple, extended, round 


Bolo 


12 





16 


(0.02) 





37 


(0.05) 


0.69 


(0.07) 


0.97 


(0.10) 


2.3 


(0.2) 


8 


77 


(1 


2) 


135 (2) 


-72 (2) 


5xl0 4 


multiple, extended, round, weak 


Bolo 


13 


















0.47 


(0.02) 


1.14 


(0.05) 


29 


44 


(0 


4) 


59 (0.5) 


33 (2) 


3xl0 5 


multiple, extended, elongated 


Bolo 


14 


















0.41 


(0.02) 


1.00 


(0.04) 


24 


44 


(0 


4) 


60 (0.5) 


-43 (2) 


3xl0 5 


multiple, extended, elongated 


Bolo 


15 





25 


(0.02) 












0.25 


(0.02) 


0.60 


(0.06) 


13 


54 


(0 


8) 


60 (0.9) 


52 (13) 


1 x 10 s 


multiple, extended, round S^J 


Bolo 


16 





21 


(0.02) 












0.21 


(0.02) 


0.51 


(0.06) 


11 


53 


(1 


0) 


58 (1.0) 


-76 (16) 


9xl0 4 


multiple, extended, round 


Bolo 


17 





13 


(0.02) 





21 


(0.04) 


0.36 


(0.06) 


0.36 


(0.06) 


0.9 


(0-14) 


9 


67 


(1 


6) 


81 (1.9) 


29 (10) 


6xl0 4 


extended, round, weak 


Bolo 


18 





25 


(0.02) 





35 


(0.04) 


0.40 


(0.06) 


0.40 


(0.06) 


1.0 


(0.15) 


16 


39 


(0 


7) 


45 (0.8) 


-64 (11) 


8xl0 5 


extended, round 


Bolo 


19 





16 


(0.02) 





35 


(0.04) 


0.59 


(0.06) 


0.59 


(0.06) 


1.4 


(0.14) 


8 


86 


(1 


4) 


93 (1.5) 


58 (16) 


4xl0 4 


extended, round, weak 


Bolo 


20 





16 


(0.02) 





37 


(0.04) 


0.64 


(0.06) 


0.69 


(0.07) 


1.7 


(0.16) 


8 


75 


(1 


2) 


102 (1.7) 


-82 (4) 


6xl0 4 


multiple, extended, round, weak 


Bolo 


21 





22 


(0.02) 












0.27 


(0.03) 


0.65 


(0.08) 


12 


44 


(0 


9) 


61 (1.2) 


88 (5) 


2 x 10 s 


multiple, extended, elongated 


Bolo 


22 





41 


(0.02) 












0.41 


(0.02) 


0.99 


(0.06) 


22 


44 


(0 


5) 


54 (0.6) 


-66 (4) 


3xl0 5 


multiple, extended, elongated 


Bolo 


23 





40 


(0.02) 












0.40 


(0.02) 


0.96 


(0.06) 


20 


52 


(o 


5) 


60 (0.6) 


75 (5) 


2xl0 5 


multiple, extended, elongated 


Bolo 


21 





32 


(0.03) 





66 


(0.05) 






0.66 


(0.05) 


1.6 


(0.12) 


16 


53 


(o 


6) 


79 (0.9) 


-1 (2) 


2xl0 5 


multiple, extended, elongated 


Bolo 


25 





12 


(0.02) 





18 


(0.04) 


0.21 


(0.07) 


0.23 


(0.08) 


0.6 


(0.18) 


8 


51 


(1 


7) 


132 (4) 


-7 (2) 


2xl0 4 


multiple, extended, round,weak 


Bolo 


26 





17 


(0.02) 





37 


(0.04) 






0.53 


(0.05) 


1.3 


(0.12) 


8 


60 


(0 


9) 


91 (1.4) 


-31 (3) 


8xl0 4 


multiple, extended, round, weak 


Bolo 


27 





17 


(0.02) 





41 


(0.04) 


0.71 


(0.06) 


1.03 


(0.09) 


2.5 


(0.2) 


11 


87 


(1 


2) 


136 (1.8) 


37 (2) 


4xl0 4 


extended, round 


Bolo 


28 





15 


(0.02) 












0.19 


(0.03) 


0.47 


(0.08) 


10 


36 


(o 


9) 


46 (1.2) 


-66 (8) 


5 x 10 5 


multiple, round 


Bolo 


29 





35 


(0.03) 












0.42 


(0.04) 


1.00 


(0.09) 


19 


33 


(o 


5) 


47 (0.7) 


-53 (4) 


lxlO 6 


multiple, extended, elongated 


Bolo 


30 





29 


(0.02) 












0.29 


(0.02) 


0.69 


(0.05) 


13 


51 


(° 


6) 


56 (0.7) 


-3 (11) 


2xl0 5 


multiple, extended, elongated 


Bolo 


31 





59 


(0.03) 


1 


24 


(0.05) 


1.92 


(0.08) 


1.92 


(0.08) 


4.6 


(0.2) 


27 


59 


(0 


4) 


90 (0.6) 


-75 (1) 


3xl0 5 


multiple, extended, elongated 


Bolo 


32 





16 


(0.02) 





37 


(0.04) 


0.63 


(0.06) 


0.76 


(0.07) 


1.8 


10-17) 


10 


94 


(1 


4) 


107 (1.6) 


-68 (9) 


4xl0 4 


extended, round 


Bolo 


33 





27 


(0.02) 












0.27 


(0.02) 


0.65 


(0.05) 


12 


55 


(° 


6) 


58 (0.6) 


90 (20) 


1 x 10 s 


multiple, extended, round 


Bolo 


34 





19 


(0.03) 












0.23 


(0.04) 


0.55 


(0.08) 


10 


40 


(1 


1) 


57 (1.6) 


4(6) 


2 x 10 s 


multiple, extended, round,weak 



Table 2 — Continued 



ID 


Flux(40") 


Flux(80") 


Flux(120") 


Total Flux 


Mass (10K) 


Peak A v 


FWHM FWHM 


PA 


<n> 


Morphology 1 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Af©) 


(mag) 


(minor,") (major,") 


(°) 


(cm ) 





Bolo 


35 





17 


(0.03) 












" 


0.20 


(0.04) 


r\ a o / r\ r\c>\ 

0.48 (0.08) 


9 


38 


(1 


2) 


56 (1.8) 


-46 (6) 


2x 10 5 


multiple, round, weak 


Bolo 


36 





13 


(0.02) 





28 


(0 


04) 


0.44 


(0.06) 


0.61 


'0.08) 


1.5 (0.18) 


9 


65 


(1 


2) 


121 (2) 


72 (2) 


5xl0 4 


extended , round 


Bolo 


37 





21 


(0.03) 





38 


(0 


06) 






0.38 


(0.06) 


u.y (u.14) 


11 


36 


(0 


8) 


75 


(1.8) 


49 (2) 


2 x10 s 


multiple, extended, round, weak 


Bolo 


38 


1 


54 


(0.03) 


2 


39 


(0 


06) 


3.33 


(0.08) 


3.33 


(0.08) 


8.0 (0.2) 


87 


47 


(0 


2) 


63 


(0.2) 


56 (1) 


2x10 


multiple, extended, elongated 


Bolo 


39 





29 


(0.03) 





59 


(0 


06) 






0.59 


(0.06) 


1.4 (U.14) 


15 


42 


(0 


7) 


76 


(1.2) 


-32 (2) 


2xl0 5 


multiple, extended, round 


Bolo 


40 





96 


(0.02) 














1.40 


(0.03) 


3.39 (0.08) 


48 


54 


(0 


2) 


63 


(0.3) 


-75 (2) 


5xl0 5 


multiple , extended , round 


Bolo 


41 





24 


(0.03) 





37 


(0 


05) 






0.37 


(0.05) 


0.9 (0.12) 


13 


44 


(0 


8) 


63 


(1.2) 


57 (4) 


2xl0 5 


multiple, extended, elongated 


Bolo 


42 


1 


54 


(0.03) 














2.54 


(0.04) 


6.1 (0.10) 


76 


49 


(0 


1) 


56 


(0.2) 


28 (2) 


2xl0 6 


multiple, extended, round 


Bolo 


43 


2 


76 


(0.03) 


4 


16 


(0 


06) 






4.46 


(0.06) 


10.8 (0.14) 


150 


47 


(0 


1) 


62 


(0.1) 


59 (1) 


2xl0 6 


multiple, extended, elongated 


Bolo 


44 





39 


(0.03) 





72 


(0 


06) 






0.72 


'0.06) 


1 t / n "i a \ 

1.7 (0.14) 


19 


45 


(0 


5) 


78 


(0.9) 


-67 (2) 


3x10° 


multiple, extended, elongated 


Bolo 


45 





68 


(0.03) 


1 


35 


(0 


06) 






1.35 


'0.06) 


6.6 (U.14J 


32 


55 


(0 


4) 


71 


(0.5) 


2 (2) 


4x 10° 


multiple, extended, elongated 


Bolo 


46 


1 


12 


(0.03) 


2 


62 


(0 


06) 






2.62 


(0.06) 


6.3 (0.14) 


52 


68 


(0 


2) 


77 


(0.3) 


41 (2) 


4xl0 5 


multiple, extended, elongated 


Bolo 


47 





93 


(0.02) 


2 


08 


(0 


05) 






2.81 


(0.06) 


6.8 (0.15) 


43 


72 


(0 


3) 


85 


(0.3) 


12 (2) 


3 x10 s 


multiple, extended, elongated 


Bolo 


48 


6 


06 


(0.03) 


9 


83 


(0 


06) 






10.59 


(0.06) 


25.6 (0.15) 


370 


39 


(0 


0) 


56 


(0.0) 


-48 (1) 


1 x 10 7 


multiple, extended, round 


Bolo 


49 


1 


13 


(0.03) 


1 


98 


(0 


06) 






1.98 


(0.06) 


4.8 (0.13) 


59 


52 


(0 


2) 


58 


(0.2) 


-49 (3) 


1 x 10 


multiple, extended, elongated* 


Bolo 


50 





45 


(0.03) 





96 


(0 


05) 






1.31 


(0.07) 


3.2 (0.16) 


22 


68 


(0 


5) 


84 


(0.6) 


-39 (3) 


2x10° 


multiple , extended , elongate^ 


Bolo 


51 





52 


(0.03) 














0.62 


(0.03) 


1.49 (0.08) 


30 


41 


(0 


3) 


58 


(0.5) 


-51 (2) 


5xl0 5 


multiple , extended , round 


Bolo 


52 





31 


(0.02) 





55 


(0 


05) 






0.79 


'0.07) 


1.9 (0.16) 


19 


44 


(0 


5) 


103 


(1.3) 


89 (1) 


2xl0 5 


multiple, extended, elongated 


Bolo 


53 





46 


(0.03) 





92 


(0 


05) 






1.37 


(0.07) 


3.3 (0.17) 


24 


61 


(0 


5) 


94 


(0.7) 


13 (1) 


2xl0 5 


multiple, extended, elongated 


r>olo 


z.a 

04 


u 


A K 


aq\ 
^O.Uo ) 


U 


yu 


/n 


Uo ) 






1 OA 


U.O / ) 


Q A 1 '7^ 
O.U (U.l I ) 


zo 


4 / 




A \ 
4 ) 


81 


(0.7) 


72 (1) 


4X 1U 


multiple, extended, elongated 


Bolo 


55 





29 


(0.02) 














0.37 


(0.03) 


0.89 (0.07) 


13 


46 


(o 


6) 


59 


(0.8) 


39 (4) 


2xl0 5 


multiple, extended, elongated 


Bolo 


56 





15 


(0.02) 





32 


(° 


04) 


0.50 


(0.06) 


0.59 


(0.07) 


1.4 (0.16) 


7 


78 


(1 


4) 


104 


(1.9) 


76 (5) 


4xl0 4 


extended, round, weak 


Bolo 


57 





30 


(0.02) 





51 


(° 


04) 






0.51 


(0.04) 


1.2 (0.10) 


16 


42 


(o 


6) 


56 


(0.7) 


-80 (4) 


4xl0 5 


multiple, extended, round 


Bolo 


58 





26 


(0.02) 














0.33 


(0.03) 


0.78 (0.08) 


19 


39 


(o 


7) 


45 


(0.8) 


-47 (10) 


6xl0 5 


multiple, extended, elongated 


Bolo 


59 





30 


(0.02) 





52 


(o 


04) 






0.52 


(0.04) 


1.24 (0.09) 


18 


46 


(o 


5) 


56 


(0.6) 


50 (5) 


3xl0 5 


multiple, extended, elongated 


Bolo 


60 





16 


(0.02) 





28 


(° 


04) 


0.37 


(0.06) 


0.37 


(0.06) 


0.9 (0.14) 


9 


44 


(1 


0) 


76 


(1.7) 


-37 (3) 


1x10 s 


extended, round,weak 


Bolo 


61 





12 


(0.02) 





22 


(0 


05) 


0.37 


(0.07) 


0.37 


(0.07) 


0.9 (0.17) 


8 


58 


(1 


7) 


104 (3) 


50 (4) 


5xl0 4 


multiple, round, weak 


Bolo 


62 





27 


(0.02) 





53 


(0 


05) 


0.80 


(0.07) 


0.80 


(0.07) 


1.9 (0.17) 


14 


65 


(o 


9) 


80 (1.1) 


-38 (6) 


1 x 10 5 


multiple, extended, elongated 


Bolo 


63 





21 


(0.02) 





44 


(0 


04) 


0.67 


(0.06) 


0.67 


(0.06) 


1.6 (0.13) 


10 


82 


(1 


0) 


89 (1.1) 


-12 (12) 


6xl0 4 


extended, elongated 


Bolo 


64 





11 


(0.02) 





17 


(0 


04) 






0.20 


(0.05) 


0.5 (0.12) 


6 


37 


(1 


2) 


82 (3) 


53 (3) 


lxlO 5 


extended, round, weak 


Bolo 


65 





67 


(0.02) 














0.76 


(0.03) 


1.84 (0.07) 


37 


41 


(o 


3) 


48 (0.3) 


4(4) 


1 x 10 s 


multiple, extended, elongated 


Bolo 


66 


1 


32 


(0.02) 


1 


83 


(0 


04) 


2.31 


(0.06) 


2.31 


(0.06) 


5.6 (0.15) 


76 


41 


(o 


1) 


44 (0.1) 


-32 (4) 


4x10 s 


extended, round 


Bolo 


67 





23 


(0.03) 





45 


(0 


05) 


0.70 


(0.07) 


0.86 


(0.10) 


2.1 (0.2) 


12 


47 


(° 


7) 


125 (1.8) 


37 (1) 


1 x 10 5 


extended, round 


Bolo 


68 





15 


(0.02) 





29 


(o 


05) 






0.36 


(0.07) 


0.9 (0.16) 


8 


43 


(1 


0) 


95 (2) 


56 (2) 


9xl0 4 


extended, round,weak 



Table 2 — Continued 



ID 


Flux(40") 


Flux(80") 


Flux(120") 


Total Flux 


Mass (10K) 


Peak A v 


FWHM 


FWHM 


PA 


(n) 


Morphology 1 






(Jy) 




(Jy) 


(Jy) 


(Jy) 


(Mq) 


(mag) 


(minor,") 


(major,") 


(°) 


(cm 3 ) 




Bolo 


69 


0.17 


(0.02) 





42 (0.04) 


0.82 (0.07) 


1.28 


(0.09) 


3.1 (0.2) 


11 


109 


(1.3) 


128 (1.6) 


59 (6) 


4xl0 4 


multiple , extended , round 


Bolo 


70 


0.35 


(0.02) 





85 (0.05) 


1.59 (0.07) 


1.96 


(0.08) 


4.7 (0.2) 


17 


94 


(0.7) 


112 (0.8) 


-63 (4) 


9xl0 4 


multiple, extended, elongated 


Bolo 


71 


0.28 


(0.03) 





69 (0.05) 




1.15 


(0.07) 


2.8 (0.18) 


13 


79 


(0.8) 


102 (1.1) 


61 (3) 


9xl0 4 


multiple, extended, elongated 


Bolo 


72 


0.34 


(0.03) 





78 (0.06) 




0.93 


(0.06) 


2.2 (0.15) 


17 


59 


(0.7) 


84 (1.0) 


38 (3) 


2xl0 5 


multiple, extended, elongated 


Bolo 


73 


0.23 


(0.02) 





53 (0.04) 


0.93 (0.07) 


1.11 


(0.08) 


2.7 (0.18) 


13 


87 


(1.0) 


114 (1.3) 


-79 (3) 


6xl0 4 


multiple, extended, elongated 


Bolo 


74 


0.38 


(0.03) 








0.38 


(0.03) 


0.91 (0.07) 


18 


51 


(0.6) 


60 (0.7) 


52 (6) 


2xl0 5 


multiple, extended, elongated 


Bolo 


75 


0.38 


(0.03) 








0.38 


(0.03) 


0.91 (0.07) 


18 


51 


(0.6) 


61 (0.7) 


-54 (6) 


2xl0 5 


multiple, extended, elongated 


Bolo 


76 


0.16 


(0.02) 





32 (0.04) 




0.43 


(0.05) 


1.0 (0.13) 


8 


64 


(1.2) 


89 (1.7) 


89 (5) 


6xl0 4 


multiple, extended, round, weak 


Bolo 


77 


0.13 


(0.03) 





24 (0.05) 


0.38 (0.08) 


0.38 


(0.08) 


0.9 (0.18) 


8 


72 


(2) 


97 (3) 


-62 (8) 


4xl0 4 


multiple, round, weak 


Bolo 


78 


0.36 


(0.02) 





72 (0.04) 




0.93 


(0.05) 


2.3 (0.13) 


19 


65 


(0.6) 


77 (0.7) 


-78 (4) 


2 x10 s 


multiple, extended, round 


Bolo 


79 


0.98 


(0.03) 








1.84 


(0.04) 


4.4 (0.11) 


43 


59 


(0.2) 


64 (0.3) 


4(4) 


6xl0 5 


multiple, extended, round 


Bolo 


80 


1.36 


(0.03) 


2 


26 (0.06) 




3.08 


(0.08) 


7.4 (0.18) 


77 


57 


(0.2) 


76 (0.3) 


-62 (1) 


7xl0 5 


multiple, extended, elongated 


Bolo 


81 


1.62 


(0.03) 








2.61 


(0.05) 


6.3 (0.12) 


86 


53 


(0.2) 


59 (0.2) 


-26 (2) 


1x10 s 


multiple, extended, round 


Bolo 


82 


0.20 


(0.02) 








0.36 


(0.04) 


0.88 (0.09) 


10 


54 


(1.0) 


62 (1.2) 


-4 (11) 


1 x 10 5 


multiple, extended, round 


Bolo 


83 


0.21 


(0.02) 





47 (0.05) 




0.47 


(0.05) 


1.1 (0.11) 


12 


60 


(1.0) 


74 (1.2) 


75 (6) 


1 x 10 5 


multiple, extended, rounds 
multiple, extended, round 


Bolo 


84 


0.25 


(0.03) 








0.47 


(0.05) 


1.1 (0.12) 


13 


61 


(0.9) 


68 (1.0) 


21 (11) 


1 x 10 B 


Bolo 


85 


0.19 


(0.02) 





40 (0.05) 




0.40 


(0.05) 


1.0 (0.11) 


11 


54 


(1.0) 


77 (1.4) 


72 (4) 


1 x 10 B 


multiple, extended, round 


Bolo 


86 


0.13 


(0.02) 





27 (0.05) 


0.43 (0.07) 


0.64 


(0.10) 


1.5 (0.2) 


9 


80 


(1.7) 


132 (3) 


-39 (3) 


3xl0 4 


extended, round,weak 


Bolo 


87 


0.13 


(0.02) 





25 (0.04) 


0.41 (0.06) 


0.45 


(0.07) 


1.1 (0.16) 


7 


68 


(1.5) 


101 (2) 


70 (4) 


4xl0 4 


extended, round, weak 


Bolo 


88 


0.21 


(0.02) 





52 (0.04) 


0.97 (0.07) 


1.43 


(0.09) 


3.4 (0.2) 


11 


94 


(1.0) 


152 (1.6) 


-39 (2) 


4xl0 4 


extended, elongated 


Bolo 


89 


0.19 


(0.02) 





44 (0.04) 


0.74 (0.06) 


0.99 


(0.08) 


2.4 (0.2) 


11 


81 


(1.0) 


126 (1.6) 


-17 (2) 


5xl0 4 


extended, round 


Bolo 


90 












0.126 


(0.015) 


0.31 (0.04) 


10 


22 


(0.9) 


33 (1.2) 


-90 (8) 




extended, round 


Bolo 


91 


0.14 


(0.02) 





24 (0.04) 


0.35 (0.06) 


0.37 


(0.07) 


0.9 (0.17) 


8 


58 


(1.5) 


105 (3) 


-79 (3) 


4xl0 4 


extended, round, weak 


Bolo 


92 


0.13 


(0.02) 





24 (0.04) 




0.24 


(0.04) 


0.57 (0.09) 


7 


54 


(1.4) 


67 (1.7) 


56 (10) 


8xl0 4 


multiple , extended , round , weak 


Bolo 


93 


0.17 


(0.02) 





38 (0.04) 


0.66 (0.06) 


0.89 


(0.08) 


2.2 (0.2) 


8 


92 


(1.4) 


111 (1.7) 


33 (6) 


4xl0 4 


extended, round,weak 


Bolo 


94 


0.16 


(0.02) 





32 (0.04) 




0.32 


(0.04) 


0.77 (0.09) 


9 


57 


1.2) 


67 (1.3) 


66 (10) 


9xl0 4 


multiple , extended , round 


Bolo 


95 


0.14 


(0.02) 





28 (0.04) 


0.51 (0.06) 


0.77 


(0.08) 


1.9 (0.18) 


8 


93 


(1.4) 


116 (1.8) 


27 (6) 


3xl0 4 


extended, round, weak 


Bolo 


96 


0.22 


(0.02) 








0.46 


(0.04) 


1.1 (0.09) 


10 


59 


(0.8) 


70 (1.0) 


8 (6) 


1x10 s 


multiple, extended, round 


Bolo 


97 


0.22 


(0.02) 








0.35 


(0.03) 


0.85 (0.08) 


11 


60 


(0.9) 


66 (1.0) 


-40 (15) 


1 x 10 B 


multiple, extended, round, weak 


Bolo 


98 


0.20 


(0.02) 








0.32 


(0.04) 


0.77 (0.09) 


9 


59 


(1.0) 


64 (1.1) 


89 (16) 


9xl0 4 


multiple , extended , round , weak 


Bolo 


99 


0.24 


(0.02) 








0.43 


(0.03) 


1.04 (0.07) 


12 


53 


(0.6) 


61 (0.7) 


75 (6) 


2xl0 5 


multiple, extended, round 


Bolo 


100 


0.41 


(0.02) 








0.81 


(0.03) 


1.95 (0.08) 


18 


59 


(0.4) 


67 (0.5) 


-7 (5) 


2xl0 5 


multiple, extended, round 


Bolo 


101 


0.21 


(0.03) 








0.35 


(0.04) 


0.8 (0.10) 


10 


48 


(0.9) 


66 (1.3) 


49 (5) 


1 x 10 B 


multiple , extended , round , weak 


Bolo 


102 


0.63 


(0.02) 


1 


41 (0.04) 




1.41 


(0.04) 


3.4 (0.11) 


31 


62 


(0.3) 


76 (0.4) 


7(2) 


3xl0 5 


multiple, extended, round 



Table 2 — Continued 



ID 


r lux^u f 


r lux^ou ) 


r lux^izu J 


Total Flux 


iviabb ^lurvj 


Pnci L- J,, 

redK s\V 


FWHM 


r vv nivi 


PA 


W 


i 

JMorpnology 






\ J y) 


W) 


(•>y) 


(Jy) 


( A/r~ \ 
\ M Q) 


(mag) 


(minor,") 


(major,") 


\ ) 


(cm - ' ) 




Bolo 


103 


1.28 (0.03) 


2.07 (0.05) 




2.42 


(0.06) 


5.8 (0.15) 


70 


46 


(0.2) 


57 (0.2) 


49 (1) 


2xl0 6 


multiple, extended, elongated 


Bolo 


104 


1.05 (0.02) 






1.42 


(0.03) 


3.43 (0.08) 


55 


46 


(0.2) 


55 (0.2) 


77 (2) 


lxlO 6 


multiple , extended , round 


Bolo 


105 


0.40 (0.02) 






0.61 


(0.03) 


1.47 (0.08) 


20 


51 


(0.4) 


64 (0.5) 


-7 (3) 


2xl0 5 


multiple, extended, elongated 


Bolo 


106 








0.50 


(0.02) 


1.21 (0.04) 


30 


52 


(0.4) 


55 (0.4) 


-87 (10) 


3xl0 5 


multiple , extended , round 


Bolo 


107 








0.48 


(0.02) 


1.17 (0.04) 


27 


53 


(0.3) 


60 (0.4) 


-76 (4) 


2xl0 5 


multiple , extended , round 


Bolo 


108 


0.18 (0.02) 






0.32 


(0.04) 


0.78 (0.09) 


10 


53 


(0.9) 


72 (1.2) 


-59 (5) 


9xl0 4 


multiple, extended, round, weak 


Bolo 


109 


0.16 (0.02) 


0.33 (0.05) 




0.39 


(0.05) 


0.9 (0.13) 


9 


50 


(1.1) 


80 (1.8) 


61 (4) 


1x10 s 


multiple, extended, round, weak 


Bolo 


110 


0.44 (0.03) 






0.44 


(0.03) 


1.05 (0.06) 


18 


53 


(0.5) 


59 (0.5) 


89 (7) 


2xl0 5 


multiple, extended, elongated 


Bolo 


111 


0.22 (0.02) 


0.47 (0.04) 


0.71 (0.07) 


0.83 


(0.08) 


2.0 (0.19) 


12 


71 


(1.0) 


101 (1.5) 


-45 (3) 


8xl0 4 


multiple, extended, round 


Bolo 


112 


0.17 (0.02) 


0.33 (0.04) 




0.45 


(0.06) 


1.1 (0.14) 


9 


60 


(1.3) 


65 (1.3) 


-70 (20) 


1 x 10 5 


multiple, extended, round, weak 


Bolo 


113 


0.18 (0.02) 


0.35 (0.05) 


0.56 (0.07) 


0.74 


(0.09) 


1.8 (0.2) 


9 


75 


(1.3) 


112 (2) 


9 (3) 


5xl0 4 


multiple, extended, round, weak 


Bolo 


114 


0.18 (0.02) 


0.38 (0.04) 




0.55 


(0.05) 


1.3 (0.13) 


10 


59 


(1.0) 


88 (1.4) 


-53 (3) 


9xl0 4 


multiple, extended, elongated 


Bolo 


115 


0.22 (0.02) 


0.37 (0.05) 


0.53 (0.07) 


0.54 


(0.08) 


1.3 (0.19) 


12 


45 


(0.8) 


96 (1.8) 


-39 (2) 


1x10 s 


extended, elongated 


Bolo 


116 


0.39 (0.02) 


0.51 (0.05) 




0.51 


(0.05) 


1.2 (0.12) 


20 


35 


(0.5) 


47 (0.6) 


46 (4) 


1x10 s 


multiple, extended, round 


Bolo 


117 


0.22 (0.02) 


0.40 (0.05) 




0.40 


(0.05) 


1.0 (0.11) 


11 


44 


(0.8) 


70 (1.3) 


76 (3) 


2xl0 5 


multiple, extended, round, \^&k 


Bolo 


118 


0.14 (0.02) 


0.29 (0.05) 




0.34 


(0.05) 


0.8 (0.12) 


8 


47 


(1.1) 


88 (2) 


4(3) 


9xl0 4 


multiple , extended , round , weak 


Bolo 


119 


0.24 (0.02) 


0.53 (0.04) 


0.91 (0.06) 


1.29 


(0.09) 


3.1 (0.2) 


13 


92 


(1.0) 


123 (1.4) 


-89 (3) 


5xl0 4 


extended, elongated 


Bolo 


120 








0.08 


(0.01) 


0.19 (0.02) 


6 


24 


(0.9) 


25 (1.0) 


-3 (90) 




multiple, extended, round, weak 


Bolo 


121 


0.18 (0.02) 


0.44 (0.04) 




0.65 


(0.05) 


1.6 (0.12) 


10 


85 


(1.1) 


94 (1.2) 


-86 (11) 


5xl0 4 


multiple, extended, elongated 


Bolo 


122 


0.38 (0.02) 


0.88 (0.04) 




1.34 


(0.06) 


3.2 (0.15) 


19 


63 


(0.5) 


91 (0.7) 


-63 (2) 


2xl0 5 


multiple , extended , elongated 



Note. — Masses are calculated according to Equation 3 from the total flux density assuming a single dust temperature of To = 10K and a dust opacity at 1.1mm of 
su mm = 0.0114 cm 2 g _1 . Peak Ay is calculated from the peak flux density as in Section 4.2. FWHM and PAs are from a gaussian fit; the PA is measured in degrees east of 
north of the major axis, (n) is the mean particle density as calculated from the total mass and the deconvolved average FWHM size. Numbers in parentheses are la uncertainties. 
Uncertainties for masses are from photometry only, and do not include uncertainties from k, Tp, or d, which can be up to a factor of a few or more. Uncertainties for the FWHM 
and PA arc formal fitting errors from the elliptical gaussian fit; additional uncertainties of 10 — 15% apply to the FWHM, and ~ 5° to the PA (determined from simulations). 

lr The morphology keyword(s) given indicates whether the source is multiple (within 3' of another source), extended (major axis FW at 2cr > 1'), elongated (axis ratio at 
4cr > 1.2), round (axis ratio at 4cx < 1.2), or weak (peak flux densities less than 5 times the RMS/pixel in the unfiltered map). 



